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ABSTRACT 


Unmanned Aircraft Systems (UASs) are critical for future combat effectiveness. 
Military planners from all branches of the Department of Defense now recognize the 
value that real time intelligence and surveillance from UASs provides the battlefield 
commander. The Operations Analysis Division of the Marine Corps Combat 
Development Command is currently conducting an Overarching Unmanned Aircraft 
Systems study to determine future force requirements. Current analysis is conducted 
through the use of the Assignment Scheduling Capability for Unmanned Air Vehicles 
(ASC-U) and several specially designed heuristics. The Unmanned Aircraft System 
Scheduling Tool (UAS-ST) combines these capabilities into one model and addresses 
several issues associated with ASC-U. UAS-ST allows the user to control all aspects of 
the UAS, define a scenario, and then generates a flight schedule over a known time 
horizon based on those inputs. All missions are assigned a user defined value and the 
total schedule value is reported. The user can then quickly change a parameter of the 
UAS, re-solve the model, and see the impact their proposed change has on the overall 
value of the schedule attained. Therefore, UAS-ST is a tool for analyzing the value of 


future changes in UAS structure. 
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THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may 
not have been exercised for all cases of interest. While every effort has been made, 
within the time available, to ensure that the programs are free of computational and logic 
errors, they cannot be considered validated. Any application of these programs without 


additional verification is at the risk of the user. 
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EXECUTIVE SUMMARY 


Unmanned Aircraft Systems (UASs) procurement is vital to the United States 
Marine Corps' (USMC) combat effectiveness in the near future. UASs are used for 
collecting intelligence, surveillance, and targeting information. They accomplish these 
missions at a much lower overall risk than conventional, manned aircraft. Military 
planners from all branches of the Department of Defense (DoD) recognize the value that 
real time intelligence and surveillance from UASs provides the battlefield commander. 
In an ongoing operation, mission requests from units in theater typically far exceed the 
capacity of available UAS assets. Demand for UAS missions is increasing as the 
capability of these platforms expands. Individual branches of DoD are scrambling to 
acquire the UASs needed to support the requirements of currently deployed units in 
combat operations. As a result, interoperability and compatibility is a major concern with 


today’s current family of UASs. 


Due to operational requirements in Iraq and Afghanistan, DoD is focusing 
procurement strategy away from force transformation and future generation weapons to 
more immediate concerns. Acquisition of UASs needed to meet operational requirements 
is the agency’s highest priority. As a result, large procurement budgets exist to fill the 
supply shortage and meet future operational requirements. Currently DoD is evaluating a 
large number of alternative UAS programs. This selection process requires an objective 
analysis of each alternative. The model in this thesis allows the USMC to accomplish 
their analysis by clearly depicting the impact of changes in system capabilities on a daily 
flight schedule. Application of the model is not limited to UASs. The model is capable 
of analyzing any asset which must be scheduled in response to user demands. Therefore, 


this model has future application in a multitude of other programs. 


The current model in use for the Marine Corps Combat Development Command 
(MCCDC) Overarching Unmanned Aircraft Systems (OUAS) study is the Assignment 
Scheduling Capability for UAVs (ASC-U) model. The design agency for ASC-U is the 
U.S. Army Training and Doctrine Command Analysis Center (TRAC). ASC-U is 


XV 


designed to support the development of an effective UAS force structure. It is a 
spreadsheet-based decision support tool primarily for allocation and scheduling of assets. 
ASC-U addresses complexities in military operations and scheduling of multiple moving 
platforms. ASC-U accepts user parameters that define a scenario then seeks to provide a 
feasible schedule for available UASs. ASC-U is the first model that enables analysts to 
address this scheduling problem effectively. ASC-U combines both optimization and 


simulation to produce a tool with unique capabilities. 


This thesis develops an integer linear programming model, UAS-ST, for 
scheduling UASs. UAS-ST allows the user to define all elements, both operational and 
performance, of the UAS via an Excel spreadsheet. A schedule generator written in 
visual basic for this application then takes these elements and generates a user-defined 
number of individual schedules. This schedule generation is done for every UAV (UAV) 
that is included in the scenario. Once the predefined number of schedules is generated, 
UAS-ST creates data files for a General Algebraic Modeling System (GAMS) model. 
GAMS, through the use of CPLEX, finds a near optimal combination of individual 


schedules to produce a complete schedule for a user designated time period. 


The initial test scenario replicates a generic scenario given in the OUAS study. 
This replication provides the simplest and most direct comparison of results with ASC-U. 
The scenario time period is twenty four hours divided into ninety six intervals of fifteen 
minutes each. This scenario does not refer to a specific country, but is designed solely to 
provide the framework for determining the value of changes in UAS structure on a given 
set of mission requests. UAS force structure for the initial scenario represents a small but 
realistic composition of systems. Tier I is comprised of three separate units with four 
UAVs and one Ground Control Station (GCS) per unit. Tier If is comprised of three 
separate units with one UAV and one GCS per unit. Tier II is comprised of a single unit 


with two UAVs and two GCSs. 


UAS-ST is significantly different from ASC-U. ASC-U utilizes an Excel 
spreadsheet and pulls the data into an Access database. UAS-ST applies an Excel 
spreadsheet, which then uses Visual Basic to transfer the data to GAMS. GAMS then 


applies the CPLEX solver to quickly optimize the schedule for all mission requests. 
XV1 


Features that require supplemental heuristics in ASC-U are incorporated directly into 
UAS-ST, eliminating the need for further processing. This produces a much quicker and 
efficient analysis of alternatives. CPLEX is highly efficient in preprocessing feasible 
solutions and reduces the run time to a matter of minutes. Overall, UAS-ST provides an 


efficient update to the model presented in ASC-U. 


UASs are critical to our nation’s and the USMC's future military combat 
effectiveness. All branches of DoD recognize the need to develop an integrated network 
of UASs. To answer this need, the Marine Corps is developing the concept of the UAS 
Family of Systems (FoS). The FoS calls for a three tier structure of UASs with 
overlapping capabilities. Currently, the MCCDC Operation Analysis Division is 
conducting an OUAS study. This thesis is a direct contribution to that study, and UAS- 
ST, provides a planning tool for development of future UAS structure. UAS-ST allows 
the user to control all aspects of the UAS, define a scenario, and then generates a flight 
schedule over a known time horizon based on those inputs. All missions are assigned a 
user-defined value and the total schedule value is reported. The user can then quickly 
change an operational or performance parameter of the UAS, re-solve the model, and see 
the impact on the overall value of the schedule. UAS-ST is a tool for analyzing the value 


of future changes in UAS structure. 
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I. INTRODUCTION 


A. PURPOSE AND OVERVIEW 


Unmanned Aircraft System (UAS) procurement is vital to the United States 
Marine Corp's (USMC) combat effectiveness in the near future. UASs are used for 
collection of intelligence, surveillance, and targeting information. They accomplish these 
missions at a much lower risk than conventional, manned aircraft. Military planners from 
all branches of the Department of Defense (DoD) recognize the value that real time 
intelligence and surveillance from UASs provides the battlefield commander. In an 
ongoing operation, mission requests from units in theater typically far exceed the 
capacity of available UAS assets. Future demand for UAS missions will increase as the 
capability of these versatile platforms expands. Individual branches of DoD are 
scrambling to acquire the UASs needed to support the requirements of currently deployed 
units in combat operations. [I[HT, 2008] As a result, interoperability and compatibility is 


a major concern with today’s current family of UASs. 


To address these issues in future systems, the USMC is developing the UAS 
Family of Systems (FoS). Each system consists of an Unmanned Air Vehicle (UAV), 
Ground Control Station (GCS), Launch and Recovery Station (LRS) and a combination 
of various sensors and payload components. The components of the FoS have several 
complementary capabilities which overlap in certain mission areas. The intent is to 
create a mix of several UASs able to support various units of different sizes and levels of 


operation. 


To assess the effectiveness of a combination of UASs, quantitative models must 
be applied to provide a reasonably accurate measurement of capability and some 
guidance for efficient employment. To this end, this thesis develops a UAS planning and 
decision support tool that takes as input a planning horizon, a fleet of UAVs and their 
individual operating limits, a list of available payloads and a list of mission requests. The 
model then provides as output an operational schedule for each individual vehicle 


indicating the missions to be accomplished in a specified time horizon and the payloads 
1 


required. This schedule provides analysts the ability to quickly determine the impact of 


changes in any system parameter on the overall value of missions accomplished. 


B. BACKGROUND 


1. Problem Statement 


Due to operational requirements in Iraq and Afghanistan, DoD is focusing 
procurement strategy away from force transformation and future generation weapons to 
more immediate concerns. Acquisition of the UASs needed to meet operational 
requirements is the agency’s highest priority. [THT, 2008] As a result, large procurement 
budgets exist to fill the supply shortage and meet future operational requirements. 
Currently DoD is evaluating a large number of alternative UAS programs. The selection 
process requires an objective analysis of each alternative. This thesis develops a model 
to help the USMC with this analysis. The model in this thesis allows the USMC to 
accomplish their analysis by clearly depicting the impact of changes in system 
capabilities on a daily flight schedule. Application of the model is not limited to UASs. 
The model is capable of analyzing any asset which must be scheduled in response to user 


demands. Therefore, this model has future application in a multitude of other programs. 


In recent comments, Defense Secretary Robert Gates, addresses wasteful or 
inefficient UAS procurement programs in a speech given at the Air Force’s Air 
University at Maxwell Air Force Base: 

Because people were stuck in old ways of doing business, it's been like 


pulling teeth. While we've doubled this capability in recent months, it is 
still not good enough. [IHT, 2008] 


The fact that the Defense Secretary is unusually blunt in his criticism of current program 


development, should serve as proof of the pressure to field future UASs. 


Zs The Marine Corps Three Tier UAS Family of Systems 


USMC Three-Tier UAS Strategy 
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Figure 1. UAS Tier Structure 


The FoS divides the various UASs into three separate tiers. The tier assignments 
provide each level of the Marine Air Ground Task Force (MAGTF) an organic, 
interoperable, integrated and tailored capability that raises the situational awareness of 
the combat unit commander through a common communication network. The current 
operational conditions in Operation Enduring Freedom and Operation Iraqi Freedom 
demonstrate the importance of UAS operations and the need for commanders at all levels 
to maintain control of their respective battle space. Therefore, tiers are defined by the 
area of interest and operational level of the supporting unit, such as a company, battalion, 
or regiment. There is some operational overlap in tier capabilities; this is intentional and 


should be considered beneficial to mission accomplishment. 


The three tiers are described in [MCCDC, 2005] as follows: 
= Tier I —- Short Range UAS 


Operationally Supports Battalion, Company, and Platoon 
Performs Reconnaissance and Surveillance Missions 
Current System: Dragon Eye, Future: Raven-B 
Endurance of 1.5 hours, Combat Radius of 5 miles 
Speed: Less than 20 knots 

Capacity for Single Payload 


* Tier II — Division or Regimental UAS 


Operationally Supports Division, Marine Expeditionary Unit 
(MEU), Regiment and Battalion 


Air vehicle: Persistent, low cost, durable, low observable, easily 
transported, shipboard compatible, target acquisition capable, 
heavy fuel engine 


Current: Scan Eagle Contract 

Endurance of 15 hrs, Combat Radius of 50 miles 
Speed: 60 knots 

Capacity for Two Payloads 


" Tier WI-—UAS 


Operationally supports Marine Expeditionary Force, Marine 
Expeditionary Brigade, MEU, Division and Regiment 


Reconnaissance Surveillance and Target Acquisition, Electro- 
Optical /Infrared Imagery (EO/IR) 


Current System: Pioneer, Future Concept - Vertical Takeoff UAV 


Air vehicle (Proposed; expeditionary, sea based, Vertical Takeoff, 
large sensor payload, multi-mission capable, weaponization) 


Sensor payload (Proposed: EO/IR/laser —_ designation, 
communication relay, Signals Intelligence, Electronic Warfare) 


Endurance of 8 hours, Combat Radius of 300 miles 
Speed: 200 knots 
Capacity for three payloads 
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FoS Tiers 





Today 2020 
Tier | Tier | 4 
e Dragon Eye ¢ Joint Manpackable UAS (Raven B) 
4 \. ° 4 systems per inf battalion 
~~ —— ¢ System = 1 GCS w/3 UAVs 
Tier Il Tier Il 


e ISR services contract (Scan Eagle) 


¢ TBD endurance system 
~ \. * 9 systems per division 
f+ System = 1GCS wi3 UAVs 





Tier Ill Tier Ill 
¢ Pioneer VUAS 
¢ 2VMUs 
4 VMUs — 
: 1 system per VMU 3 systems per VMU Feo MM 


System = 1 GCS w/ 5 UAVs 





System = 2 GCS wi 4 UAVs 


Figure 2. Projected Future FoS Tiers Structure 


Current UAS Structure 


There are several specialized terms associated with a UAS that this thesis uses 


repeatedly. To avoid any confusion, the following definitions apply: 


UAV (Unmanned Aerial Vehicle): A UAV is an unpiloted aircraft. UAVs 
can be remote controlled or fly an automated route based on pre-programmed 


information. 


UAS (Unmanned Aircraft System): UAS is the current term introduced by 
DoD and accepted by the Federal Aviation Administration to replace the term 
UAV. A UAS consists of not only the unmanned aircraft, but also the data 
link system, the launch and recovery station, and all maintenance and support 


equipment. 


GCS (Ground Control Station): A GCS is a land or sea-based system that 
allows an operator to control an unmanned aircraft. A GCS may control both 
"active" and "passive" missions. An active mission requires the operator to 
monitor the UAV, while a passive mission requires no interaction from the 


operator. 


e LRS (Launch and Recovery Site): Used to control an unmanned aircraft 


during the initial and terminal phases of flight. The LRS can also function 


as a GCS when not recovering UAVs. 


e Payload: Future system requirements call for a system of “payloads,” 


each of which provides a custom capability. These payloads will be 


uniform in size and connectivity to allow for rapid configuration of 


mission specific profiles. 


e Mission Package: Consists of the combination of payloads loaded on an 


individual UAV which determines the mission capability of the UAV. 


Required payload capacity is a key element of the current study. 


4. Current UAS Missions Types 


Following is a table of terms associated with current UAS missions. These 


missions represent general categories and are derived from the Overarching Unmanned 


Aircraft System (OUAS) study [OUAS, 2007]. 


Reconnaissance, Surveillance, Target 
Acquisition (RSTA) (EO/IR) 


ISR asset for routine day and night time 
operations 





Reconnaissance, Surveillance, Target 
Acquisition (Synthetic aperture radar) 


ISR asset for dense vegetation and poor 
weather conditions 








Signals Intelligence 


Sensor(s) designed for passive collection of 
signals 








Air Vehicle Communications Link Relay 


Communications Relay 


Relays information and instructions through 
one UAV to another 


Relays voice and data between ground points 








Strike 


Weapons enabled kinetic destruction of a time 
sensitive target 











Electronic Warfare 








Active denial of radio frequency 








Table 1. | General Description of Mission Types 





C. SCOPE AND LIMITATIONS 


The intent of this thesis is to provide a model for use as an analytical tool in 
determining the future force structure for UASs. The model allows the user to quickly 
change both operational and performance parameters for each UAS tier. Once 
parameters are set, the model quickly generates a near-optimal schedule. This rapid 
schedule generation allows an analyst to see the impact of changing UAS capabilities on 


a given daily schedule. 
D. THESIS ORGANIZATION 


Chapter II provides a discussion of the Assignment Scheduling Capability for 
UAVs (ASC-U) the current UAS evaluation tool in use by the Operation Analysis 
Division (OAD), MCCDC. Chapter III describes the optimization model and the stack 
based enumeration heuristic used to solve it. Chapter IV provides a detailed analysis of 
the ASC-U results and compares it to the current model. Chapter V is devoted to 


conclusions and recommendations for future research. 
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I. ASSIGNMENT SCHEDULING CAPABILITY FOR UAS 


A. INTRODUCTION 


The current model in use for the OUAS study is the ASC-U model. The agency 
responsible for the design of ASC-U is the U.S. Army Training and Doctrine Command 
Analysis Center (TRAC). ASC-U is designed to support the development of an effective 
UAS force structure. This chapter describes the basic structure of the ASC-U model. It 
begins with a description of the inputs accepted by the model then describes its 
implementation in the study. The chapter then discusses limitations of ASC-U which are 
addressed by this thesis. The majority of the information given here is summarized from 
the OUAS report [OUAS, 2007]. This chapter also includes a review of current 


literature. 


B. BASICS OF ASC-U 


1. Model Description 


ASC-U supports the development of an effective UAS force structure. It is a 
spreadsheet-based decision support tool used primarily for allocation and scheduling. 
ASC-U addresses the complexities in military operations and scheduling of multiple 
moving platforms. ASC-U accepts user parameters that define a scenario and then seeks 
to provide a feasible schedule for available UASs. As stated in the ASC-U Analyst 
Manual: 

ASC-U provides a solution to the following problem: Given a scenario 

that specifies the number of each type of UAV, initial UAV locations, and 

UAV performance characteristics, determine the number of missions that 

can successfully be completed and the schedule for each UAV. The 

solution must consider GCS locations and capacities, remote viewing 


terminal requirements, and communication platform footprint and 
capacities. [Ahner, 2006] 


ASC-U is the first model that enables analysts to address this scheduling problem 
effectively. It combines both optimization and simulation to produce a tool with unique 


capabilities. 


2. ASC-U Implementation 


The ASC-U model works by allowing the user to design a scenario which consists 
of a set of mission requests. A mission request consists of a specific, required UAS 
capability at a specific geographic location, for a set amount of time. The location of the 
mission remains fixed once it is assigned. The allocation tool uses the UAS capability 
data and mission data it is given to create a feasible schedule which will accomplish as 
many mission requests as possible. ASC-U is deterministic. For a given input, it will 
always produce the same schedule with the same measures of performance. Specifically, 
ASC-U uses a deterministic algorithm to optimize over a given finite time horizon to 


obtain near-optimal UAS mission area assignments. 


ASC-U allows the user to define several different input parameters. Mission 
requirements are the most important parameter. Mission requirements consist of a 
coordinate location, payload requirement, length of mission and mission priority. 
Mission priority functions as a selection criterion. A user may establish mission 
precedence by assigning a higher value to a specific mission type. UAS parameters that 
may be entered by the user are payload details, GCS attributes, and UAV attributes. GCS 
attributes consist of coordinate location, control limits, and unit assignment. UAV 
attributes include speed, operating time, combat radius, launch site, and total time 


available. 
3s Model Capabilities and Limitations 


ASC-U is the first attempt to create a specific tool that involves all aspects of the 
UAS family of systems. As part of its support program, TRAC-Monterey publishes an 
ASC-U user’s manual. For its use in the OUAS study, some settings are different than 
the manual’s recommended settings. Different settings are required due to the small size 


of the original test scenario chosen. ASC-U is designed to model several thousands of 
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missions over a long time period. The scenarios used in the OUAS study are much 


shorter in duration and involve only a few hundred missions. 


The mission requirements are the most important input data because they have the 
most direct effect on the schedule generated. ASC-U attempts to complete as many 
mission hours as possible. The objective function has a significant drawback. It leads to 
preferential assignment of missions that are close to the launch and recovery site because 
they allow for more follow-on missions to be accomplished. The model schedules as 
many close-in missions as possible and foregoes farther outlying missions. Therefore, 
mission priority, a specific input, allows the user to set a precedence level. However, in 
some instances the user is forced to artificially inflate the value of a mission to ensure it 


is scheduled. 


All inputs are usually entered in Excel spreadsheets, and then read into Access. 
ASC-U can take input from either Excel or Access. The Excel inputs are placed in an 
Access data base as the model begins its run. The output is also stored in an Access data 
base. Several tables are produced in the output. The most important table for the OUAS 
study is mission coverage. Mission coverage is broken down into requesting unit, 
mission type and UAV type. In ASC-U every possible payload combination is 
enumerated as a mission package. Mission packages are then assigned to UAVs. 


Mission package usage is recorded; therefore individual payload usage is not available. 


The objective of the first OUAS study is to provide as many hours of UAS 
support as possible. For ASC-U, the two key factors are the optimization interval and the 
time horizon for scheduling the UAV. The optimization interval controls how often 
ASC-U runs its routine to assign missions. The time horizon controls how far forward in 
time ASC-U will look to assign the UAV to a mission. The optimization interval is 
crucial because if the interval is set too long then the UAS will remain idle instead of 
performing another mission. If the interval is set too low the model has difficulties and 
exhibits odd behavior. An interval of 12 minutes is best for the scenarios used in the 
OUAS study [OUAS, 2007]. Time horizon is critical because ASC-U will only schedule 
a UAV once during the given time period. Therefore, if the horizon is set too long, ASC- 


U can see a high value mission at the end of the period and keep a UAV idle the entire 
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period waiting for that mission. A UAV can complete another mission and still complete 
its priority mission. On the other hand, if the time horizon is set too low then ASC-U 
misses high priority missions because it cannot see far enough out. The final horizon for 


the OUAS study is 6 hours for Tier II and 4 hours for Tier III. 


Because of these shortcomings, several workarounds and new heuristics are used 
to enable ASC-U to find schedules that exhibit required behaviors. The SUPER 
MISSION ability allows a UAV to fulfill multiple missions if they are within its operating 
range. Three optional heuristics are also used in the OUAS study. They are implemented 


in this order: 


EARLY RETURN (Go Home vs. Stay) - When a UAV completes an assigned 
mission, a value of TRUE allows the UAV to return to base if it can do so and still make 
it back on station in time for its next mission. A value of FALSE forces the UAV to fly 


until its operating limit is met. For the OUAS study, this is set to FALSE. 


SECONDARY AREAS (Go Get More Value From Another Mission) - If a UAV is 
scheduled for more than one mission in the same location with a time gap in between 
them, a value of TRUE allows the UAV to perform another mission in between as long as 
it is available at the start of its previously scheduled mission. A value of FALSE will 
force it to remain on station until the start of its next mission. For the OUAS study, this 


was set to TRUE. 


APPENDED AREAS (Done, Is There Another Mission?) - If a UAV has 
completed its assigned missions, a value of TRUE allows it to use its remaining time to 
find another mission. A value of FALSE will not allow additional missions. For the 


OUAS study, this was set to TRUE. 


C; LITERATURE REVIEW 


Many recent studies attempt to shape some aspect of the future design of UAS 
force structure. An unintended consequence is that most of these studies focus only on 
some specific technical portion of the overall problem. The Deputy Commandant, 


Aviation and the Deputy Commandant Combat Development and Integration are 
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sponsoring the OUAS study. Their intent is to analyze the future USMC UAS force 
structure to determine how to best meet the needs of the MAGTF. The initial phase of 
OUAS is being conducted by OAD MCCDC. The model given in this thesis is intended 
for use in the follow on analysis of OUAS. Therefore, the initial report by OAD is 


critical to this thesis, as it largely determined the requirements for the model. 


Tutton [2003] deals with the optimal placement of a unit’s sensing assets. He 
presents a methodology for finding the most beneficial mix and allocation strategy for an 
individual unit’s sensors for a given threat scenario. Doll [2004] takes the model 
developed by Tutton [2003] and translates it into a programming language for easier 
simulation. She refines many of the constraints in the original model to make for a much 
more realistic simulation. Finally, Zacherl [2006] deals specifically with reactive aircraft 
scheduling. The thesis develops a model which reviews a current air tasking order and 


then rapidly reassigns aircraft to new targets as they become available. 


Finally a large number of commercial information sources address current topics 
in UAS development and model optimization. When possible, these sources are used for 


the sake of currency or accuracy. 
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HW. OPTIMIZATION MODEL FOR THE SCHEDULING OF 
UNMANNED AIRCRAFT SYSTEMS 


A. INTRODUCTION 


This chapter develops the mathematical programming model, Unmanned Aircraft 
System Scheduling Tool (UAS-ST). UAS-ST applies an optimization based approach to 
analyze the best mix of future UAS mission capabilities. The main goal of this thesis is 


to provide a model, which will be applied to the current OUAS study. 


The model allows the user to define all elements of the UAS via Excel 
spreadsheet. The schedule generator then takes these elements and generates a user 
defined number of individual schedules for every UAV included in the scenario. This 
schedule generation is done via Excel and Visual Basic. Once the predefined number of 
schedules is generated, the data is read into a General Algebraic Modeling System 


(GAMS) program to find a near optimal combination of individual schedules. 


B. AN INTEGER PROGRAM TO OPTIMIZE UAS SCHEDULING 


The following integer linear program (ILP), UAS-ST, attempts to find the 
consolidated UAS schedule with the highest overall total value. First the model is 
presented, then the detailed input required to optimize the objective function is discussed. 


Once this is complete, instances of UAS-ST are generated to demonstrate its function. 


1. Indices 

veV Tier levels [3] 

heH Mission types [2] (H= {passive, active }) 
geEG Ground Control Station (GCS) [~10] 
leL Launch and Recovery Site (LRS) [~6] 
meM Missions [~150] (alias m’) 
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peP 
teT 
ses 


m'eP 


m 


(m,m')eE 


(s,m,t) eA 
(s,v.t)eB 
(s,1,t)eK 
(s,c,t)€Q 
(s,1) SL 
(m,g)<eMG 
(m,h) « MH 


(g.v)e GV 


Tn 
num_UAVS, ¢ 


num_payet 


Payload module types [~10] 
Time periods [96] 
UAV employment schedules 


Mission prerequisites: mission m cannot be covered in any time 


period ¢ unless at least one mission m’ in P,, is also covered in time 


period t 


Pairwise exclusive missions: mission m cannot be covered in any 


time period that mission m’ is covered. EGCMxM. 
schedule s covers mission m in time period tf 
schedule s uses UAV in tier level v in time period tf 
schedule s uses LRS / in time period t 

schedule s carries payload of class c in time period t 
schedule s uses LRS / 

mission m can be covered by GCS g 

mission m is of mission type h 


GCS g can support a UAS in tier v 


Given Data [units] 


Value per time period of mission m (=total value of m divided by 


length of m) 
Number of time periods of required coverage for mission m 
Number of UAVs in tier level v available in time period ¢ 


Number of payloads of class c available in time period tf 
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BCs _ cap, Capacity (in UAVs) of GCS g in time period ¢ for missions of 


type h 


Irs_capit Capacity (in UAVs) of LRS / in time period t 


lengthn Maximum number of periods mission m can be covered 
Kids =1 if schedule s uses LRS / in time period t 
3s Decision Variables 
X, =1 if schedule s flown by some UAS [binary] 
oat =1 if mission m covered by a tier v UAS in period ¢ 
eee =1 if mission m supported by GCS g for a tier v UAS in time 
period ¢ [binary] 
a Total dwell time on mission m [time periods] 
LOAD, , ; Payloads of class c for tier v sited at LRS / [cardinality] 
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4. 


Formulation 


max > val,,D,, 
m 


S.t. 


> XxX, <num_uavs,,, 


si(s,v,t)eB 
(s,J)eSL 


» X,<LOAD,,, 


ek 


>) LOAD <num _ pay. 
I 


h 
DS Wog35 < gcs_cap,, 


m,v:(m,g )EMG 
m,h)eMH 


g.v)eGV 





X , <Irs_cap,, 
si(s,Lt)eK 


Gare Ss: x 


s:(s,m,t)eA 
(s,v,t)eB 


DY ue 


a s >; Wrrsices 


g:(m,g)eMG 
(g,v)eGV 


Deore + ae = 1 
D. 2f, 

0<D,, <length,, 
X,€ {0, 1} 

a= {0,1} 


Wriver € {051} 


m,V,g,t 


Vi,v,t 
Ve,v,l,t 


Vov 


Vg,t,h 


Vm: P, #0,Vv,Vt 


Vm: P, =O,Vv,Vt 


V(m,m') e€E,Vt 
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(RO) 


(RI) 


(R2a) 


(R25) 


(R3) 


(R4) 


(RS) 


(R6) 
(R7a) 


(R7b) 


(R8) 


(RY) 

(R10) 
(R11) 
(R12) 
(R13) 


5. Discussion 


The objective function (RO) calculates the value of all mission time periods 
covered. Constraints (R1) limit the number of active UAVs in each tier, in each time 
period, based at each LRS, by the number of UAVS available. Constraints (R2a) limit the 
number of a specific payload class an tier flown from an LRS in each period to the 
number of payloads of that class and tier assigned to that LRS, and constraints (R2b) 
limit the total number of payloads of each class and each tier assigned over all LRSs by 
the total number available in that class. Constraints (R3) limit the number of missions, 
by type supported by a GCS in each time period. Constraints (R4) limit the total number 
of UAVs launching or recovering at a LRS by the capacity of the LRS in a given time 
period. Constraints (R5) controls whether or not a mission is accomplished in any time 
period. Constraints (R6) limit the total time spent on a mission to the number of time 
periods covered. Constraints (R7a) address prerequisite missions such as AV_CLR 
required for long range Tier III missions, and (R7b) address line-of-sight issues from 
each GCS for missions that do not have other prerequisites. Constraints (R8) prevent 
scheduling of mutually exclusive missions by UAS of any tier. Constraints (R9) force 
“required” missions to be covered for the number of time periods required. Constraints 
(R10) limit the total dwell time on a mission to between zero and the total amount of time 
requested for the mission. Constraints (R11-R13) restrict schedule assignment, mission 


coverage, and GCS assignment to be binary decisions. 
6. Route Enumeration and Data Development 


UAS-ST requires a significant amount of data processing to produce an optimal 
schedule; most of this effort is performed by a stack-based enumeration routine that 
generates all feasible schedules for each tier and LRS in the given scenario. The schedule 
generator, or simply generator, requires two primary data structures for computation of 
feasible schedules: a stack, PATH(), containing a current feasible list of missions to be 
accomplished, and an array, ON_STACK(), that indicates whether each mission is 


currently on the stack of missions. Various data tables are pre-computed to aid in 
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determining feasibility, such as flight times between airfields and missions, and between 


pairs of missions, payload requirements by mission, etc. 


For each tier, v, and LRS, /, such that there is at least one UAS of tier v located at 
J, the generator starts with an empty mission list PATH(), (representing the trivial action 
of launching, and then immediately recovering, a UAS of tier v from J), and builds 
feasible mission lists exhaustively by adding one mission at a time to the end of the 
current mission list represented in PATH(). A mission is only added to the end of 
PATHO if it is within range of a GCS capable of controlling the appropriate tier, if the 
first period of the mission that can be accomplished (based on the current mission list) 
allows the UAS to return back to its LRS before running out of operational time 
available, and if the required list of payloads to accomplish all the missions in PATH() 
can fit on one UAS of tier v. If a mission does not meet these requirements, it is 
discarded and the next mission from the overall mission list is considered. Other rules for 
calculating feasible extensions to a current feasible PATH() can be incorporated easily, 
but these three capture the primary constraints on feasibility. For example, we also check 
to see if a mission is already on the stack PATH(), and discard those to prevent missions 
being revisited. However, if the capability to revisit missions in the same schedule is 
required, this test can be removed. Of course, more feasible schedules will be generated 


if this is done. 


Once a mission is added to PATH), all missions are considered again for the next 
empty slot at the end of PATH(). In this manner the generator provides a depth-first 
exploration of all feasible mission lists. If all missions have been ruled out for a slot, the 
stack is “popped,” the previous slot is then considered again, and the next mission in the 
list takes the current mission’s place; this replacement is repeated until a feasible mission 
for that slot is found. If no more feasible missions are found for that slot, we “pop” the 
stack again. When the first slot is finished, we have enumerated all missions from / using 


tier v assets. 


Every feasible set of missions in PATH() is recorded by incrementing the number 


of feasible paths found, s, and then adding the appropriate elements to the sets 
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(s,m,t)eA, (s,v.t)eB, (s,c,t)€Q, and (s,/)eSLto define schedule s. See 
Appendix B for the details of the algorithm. 


We set an upper bound on the total number of schedules that can be generated, so 
as not to create unsolvable ILPs. We also provide a limit on the number of missions 
enumerated per schedule, and the generator implicitly calculates a limit on the number of 
schedules to generate for each tier-LRS pair, to allow for the generation of a nonempty 
list for each such pair regardless of the order in which schedules are generated for those 
pairs. So, for example, we might set a limit of 120,000 total schedules, and four missions 
per schedule. The generator calculates how many tier-LRS pairs are feasible, and if, say, 
there are six such pairs in the scenario, then the generator will generate no more than 
20,000 schedules for each pair considered, calculated cumulatively. (Specifically, in this 
example if the first pair does not generate all 20,000 of its allotted schedules, then the 
generator can generate up to 40,000 total schedules for the first two pairs combined.) 
When these enumeration limits are reached for a tier-LRS pair, the stack is cleared out 


and the generator moves to the next tier-LRS pair. 


This procedure can generate tens- or hundreds-of-thousands of feasible schedules. 
Each such schedule is associated with a decision variable, X,, in the ILP. The limits on 
enumeration therefore restrict the schedules available, and, consequently, the resulting 
ILP will be a restriction of the model that considers every feasible schedule. The more 
schedules that are generated (through increasing the limit on missions per schedule, or the 


total schedules generated, etc.) the better the final schedule found by the integer program. 
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IV. COMPUTATIONAL RESULTS 


The OUAS study generic unclassified scenario is replicated with the intent to test 
the results provided by UAS-ST [OUAS, 2007]. All computations are performed on a 
Pentium 4 desktop computer at the Naval Postgraduate School with the use of the GAMS 
solver CPLEX [2008]. 


A. INITIAL TEST SCENARIO 


The initial test scenario replicates a generic scenario given in the OUAS study 
[OUAS, 2007]. This replication provides the simplest and most direct comparison of 
results with ASC-U. The scenario time period is twenty four hours divided into ninety 
six intervals of fifteen minutes each. This scenario does not refer to a specific country, 
but is designed solely to provide the framework for determining the value of changes in 
UAS structure on a given set of mission requests. UAS force structure for the initial 
scenario represents a small but realistic composition of systems. Tier I is comprised of 
three separate units with four UAVs and one Ground Control Station (GCS) per unit. 
Tier II is comprised of three separate units with one UAV and one GCS per unit. Tier III 


is comprised of a single unit with two UAVs and two GCSs. 


Table 2 provides a summary of the given UAS force structure for the generic scenario. 
































Number of: | Tier I | Tier II | Tier I 
Systems IZ B 2 
UAVs per system 3 3 4 
GCSs per system i 1 2 
UAVs controlled byGCS_| 1 2 2 





Table 2. Number of UAS Assets Available in Generic Scenario 


The initial test scenario consists of placing units from different UAS Tiers at 
specific ranges to test the scheduling constraints of the model. The tier III squadron is 


centered at an airfield in support of a regiment. The squadron consists of two UAVs and 
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two GCSs at this location. These are the only tier II assets available for our test 
scenario. The mission areas for the tier III UAS are obtained by plotting points in the 
cardinal directions at or beyond one hundred and fifty nautical miles, (beyond tier II 


range). 


The three tier II units are deployed in a triangular pattern around the tier III 
location in direct support of individual battalions. The tier II units are separated from 
each other by a distance of at least 20 nautical miles. As before, specific mission areas 


are chosen for each tier II unit. 


Tier I units are collocated in the same manner as tier II, but are in direct support at 
the company level. Each tier I unit consists of four UAVs and a single GCS. Mission 
areas are selected with a special emphasis on their range. Once the areas are designated, 
they are used repetitively to generate multiple missions. The use of repetitive mission 


areas makes it much simpler to detect mission assignments beyond the range of a specific 


tier. Tier I and II missions are within range of tier II assets. 





Figure 3. Plot of OUAS Scenario Taken from OUAS 
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B. UAS-ST INITIAL RESULTS AND ANALYSIS 


Data for the initial scenario is entered via a mission request worksheet in Excel. 


The user enters mission requests for UAS coverage in the following manner: 


Mission End Start |Finish| Periods [Per Time 
Class Time|Requested|Period/Period|Req i 
m1 


N323934|W1144850/AV_CLR}0:00/8:00]| 8:00 1 32 32 5 




































m2 N322947IW1144606] ISR {0:30/2:30] 2:00 3 10 8 10 
m3 N322825]W1144547| ISR {1:00/4:00] 3:00 5 16 12 10 
m4 N322648IW1144730| ISR |1:30/2:30 1:00 7 10 4 10 
m5 N322549IW1144914) ISR {2:00/4:00]| 2:00 ) 16 8 10 
m6 N324012IW1142426| ISR {2:30/5:30] 3:00 11 22 12 10 
m7 N3240211W1142011] ISR {3:00/4:00 1:00 13 16 4 10 
ms N324055]W1142318) ISR {3:30/5:30] 2:00 15 22 8 10 
m9 N323922IW1142442) ISR {4:00/7:00| 3:00 17 28 12 10 
mi0 _|N323848)W1142056} ISR |4:30/5:30 1:00 19 22 4 10 





Table 3. Input Worksheet for Mission Requests 


Waypoints for the mission areas are placed at specific ranges to test the 
scheduling constraints of the model. These waypoints are then used to generate a table of 
one hundred and fifty mission requests. The requests are broken down by mission areas. 
Mission requests one through forty focus on the tier I mission areas. Requests forty one 
through one hundred and three focuses on tier II and the remaining requests are tier III 
mission areas. See Appendix A for the complete scenario worksheet used in analysis of 
UAS-ST. The model constraints allow a higher tier vehicle to perform a lower tier 
mission if that is the highest value mission available at the time and the vehicle is within 
range. However, the longer range of each subsequent tier’s mission area make it 
infeasible to schedule a lower tier UAS. The sequential mission numbers enable the user 


to quickly identify any infeasible mission assignments. 


Mission classes are defined according to the mission categories given in the 
OUAS study. Values for the individual missions are chosen arbitrarily. Selection of the 
assigned mission values is critical because they have the ability to skew the results by 
inflating the value of a specific category. In the scenario, active missions are given 


greater priority than passive. Of the active subset, missions which are unique to a 
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specific tier are given the highest priority. Missions which are redundant to all tiers are 
given the lowest priority. Table 4 outlines the missions which each tier can perform and 


their associated value in the scenario. 


Tier | Tier Il Tier Ill 


ISR 
AV_CLR [5 
CR 
LSR_P 
SIGINT 
EW 
STK 
LSR_D 
Table 4. Assigned Mission Values by Tier 


























The requested mission time is defined through the start and end times, which are 
entered by the user. UAS-ST then converts that time into the appropriate number of 
fifteen minute time periods. Once a UAS is assigned in a given time period it may not be 
reassigned until the next time period. Times are started at the beginning of the first time 


period and continue to the end of the final period. 


Another element which greatly affects the optimal value achieved is the number 
of mission payloads available. See Figure 4 for a sample of the payload worksheet. By 
manipulating the number of payloads of a given type available to a UAS tier, it is 
possible to test the scheduling constraints and see very quickly if payload availability was 
a limiting factor. This is highly beneficial in determining future requirements for UAS 


structure. 
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ID Payload Type Class Available | Tier 
pl EO/IR ISR 9 
EO/IR-targeting ISR_T 
EO/IR-targeting ISR_T 
SAR w/MTI ISR_SAR 
SIGINT SIGINT 
SIGINT SIGINT 
Mine Detection ISR_MD 
Pointer LSR_P 
Pointer LSR_P 
Rangefinder LSR_RF 
Designator LSR_D 
AV_CLR AV_CLR 
AV_CLR AV_CLR 
AV_CLR 









































CR 
CR 
Strike 
EW 








Nh pm | ms fos | Oo] oO Fm ms Ss Ss st Ss Sst st Ii 








ine) 


Figure 4. Example of Payload Worksheet 


System capabilities and limitations are drawn directly from the OUAS study. The 
study seeks to define UAS tier requirements. Therefore, those assumptions are critical to 
any analysis that is conducted. The tier parameters as given in Figure 5 are used in all 
testing. Operational time limit and maintenance time are not known at this time. Those 


columns are for future expansion when the data becomes available. 


Max Cruise Operational | Payload Operational |Maintenance 
UAS Tier | Endurance (hrs) | Speed(knots) | Radius(nm) | Capacity [Time Limit(hrs) Time 


T1 1.5 17 5 1 
T2 15 60 50 2 
T3 8 200 300 3 


* All Values Taken from Table ES-2 of Overarching UAS Study 21Nov2007 





Figure 5. OUAS Study Future Tier Performance Parameters 


Once all data for the scenario is entered via Excel spreadsheet, UAS-ST is run 
from the graphical interface given in Figure 6. The interface, known as the dashboard, 
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allows the user to specify several key parameters. First the user specifies the start time, 
the time periods in horizon, and the number of periods per hour. For the mission 
generation portion, a user specifies the maximum number of schedules that UAS-ST can 
generate and the maximum number of missions per schedule that a single UAV can 


perform. These two parameters greatly affect the outcome of the model. 


As the number of maximum missions per schedule increases, the number of 
possible schedules rapidly increases as well. Setting the number of maximum schedules 
too low prevents a full enumeration of all possible schedules and is a constraint on the 
model. Setting the number of maximum schedules too high allows a full enumeration of 
all possible schedules and results in a longer processing time. See Figure 6 for an 
example. All tested combinations satisfy the requirement for a quick total solution time. 
Once all parameters are set, selecting the BUILD function initiates VBA to build the 
designated number of flight schedules. After the build is complete, selecting SOLVE 
initiates GAMS to read in the data, apply the selected solver, and generate a near optimal 
solution. When GAMS reaches a solution, selection of RESULTS displays the solution in 


the format seen in Figure 7. 


UAS Planner v1.01 2-Jun-08 
Start Time 0:00 Build | 
Time Periods in Horizon 96 Bone | 
First Planning Period 1 
Last Planning Period 96 Results | 
Periods per hour 4 
Max Schedules 75000 
Max Missions per Schedule 3 Routes: 67324 
Solver CPLEX Status: 
OPTCR 0.05 
RESLIM 3600 

Figure 6. Sample of User Interface (Dashboard) 
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Total Total 


Start End Time Periods | Time Periods | Per Time 
Time Time Requested Assigned Period |Achieved 23340 












0:00 8:00 0 10 
m2 0:30 2:30 0 5 
m3 1:00 4:00 0 5 
m4 1:30 2:30 0 5 
m5 2:00 4:00 0 5 
m6 2:30 5:30 6 5 
m7 3:00 4:00 1 5 
ms 3:30 5:30 4 5 
m9 4:00 7:00 7 5 
m10 4:30 5:30 0 5 
m11 5:00 7:00 5 5 
mi2 5:30 8:30 9 5 
m13 6:00 7:00 0 5 
m14 6:30 8:30 0 5 
m15 7:00 10:00 6 5 
m16 7:30 8:30 0 5 
m17 8:00 16:00 0 10 
m18 8:30 11:30 0 5 
m19 9:00 10:00 0 5 
m20 9:30 11:30 0 5 
m21 10:00 13:00 12 12 5 

Figure 7. Partial Display of Optimal UAS-ST Schedule 


Figure 7 demonstrates the format in which the solution given by UAS-ST is 
displayed. All missions are sorted sequentially by their assigned mission number. The 
total value achieved for each mission is given in the right column. The total value 
obtained by the entire schedule is given in the top right corner. This format enables the 
easiest comparison between runs because it quickly lets the user see which missions are 
scheduled and for what portion of the requested time. If certain missions are not being 


scheduled, the trend is easy to detect. 
C. COMPARING UAS-ST TO ASC-U 


ASC-U is the model originally used in the OUAS study. It is highly useful in the 
analysis of UAS future requirements with several strong attributes. However, there are 
some areas which require modification to meet the needs of the OUAS study. The intent 
of this thesis is to develop a user friendly model that addresses those shortcomings and 


meets the requirements of OAD. 
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For ease-of-use the individual elements of the UAS are separated into different 
worksheets. A separate spreadsheet for the UAS tier- specific performance requirements 
is also included. This spreadsheet allows the user to verify the basic parameters of all 
UAS tiers in one chart. The user can easily modify a parameter and re-run the model to 
analyze the effect. The model applies these values to the system constraints when 


optimizing the schedule. 


UAS-ST takes a different approach from ASC-U when addressing the 
optimization of payloads. In ASC-U all possible combination of payloads are 
enumerated and then designated as mission packages. UAS-ST treats the payloads as 
individual units, which allows for tracking of specific payload utilization. OAD indicated 
this was necessary for determining future requirements. The payload data is further 
separated by creating a mission class worksheet and a payload worksheet. The mission 
class worksheet defines all current missions for the UAS. It allows the user to designate 
a mission as either active or passive and which tier is capable of covering that mission. 


The payload worksheet allows the user to designate the quantity available. 


A major difference between the two models is in the handling of the mission data. 
In UAS-ST all information is consolidated, both input and output, in the same worksheet. 
The user builds the missions page in the same manner as a flight schedule and assigns a 
specific number to an individual mission. UAS-ST sorts all output based on mission 
number, which makes for easy tracking of a specific mission. ASC-U uses either the 
mission number or the coordinates of the mission location to sort output. This makes 
analysis of results extremely time consuming and requires additional sorting. ASC-U 
requires the use of a heuristic, SUPER MISSIONS, to designate a mission as mandatory 
by assigning it extra value. UAS-ST addresses this in the mission input page by allowing 
the user to designate any portion of a mission request as required. Setting this 


requirement serves as an additional constraint on the optimization model. 


ASC-U utilizes an Excel spreadsheet and pulls the data into an Access database. 
UAS-ST applies an Excel spreadsheet, which uses Visual Basic to transfer the data to 
GAMS. GAMS then applies the CPLEX solver to quickly optimize the schedule for the 


mission requests it is given. The required features addressed by the supplemental 
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heuristics in ASC-U are incorporated directly into UAS-ST to eliminate the need for 
further processing. This allows for a much quicker and efficient analysis of alternatives. 
CPLEX is highly efficient in preprocessing the feasible solutions and reduces the run 


time to a matter of minutes. Overall, UAS-ST provides an efficient update to ASC-U. 
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V. CONCLUSIONS AND NEW OPPORTUNITIES 


A. SUMMARY 


This thesis develops a scheduling tool, UAS-ST as a replacement for ASC-U. 
UAS-ST addresses several areas which require additional processing in ASC-U. UAS- 
ST allows a user to generate a daily flight based on a given set of operational and 
performance parameters. The true strength of UAS-ST is that a user can change a single 
parameter or some combination of parameters and quickly rerun the model to see the 
effect on a flight schedule. UASs are critical to our nation’s future combat effectiveness. 
As a result, all branches of DoD recognize the need to develop an integrated network of 
UASs. MCCDC OAD is currently conducting an OUAS study with the goal of defining 
future UAS requirements. This thesis is a direct contribution to that study, and the model 


developed herein, provides a planning tool for the development of future UAS structure. 
B. OPERATIONAL INTRODUCTION 


UAS-ST is currently intended for use in the second phase of the OUAS study. It 
will replace ASC-U and provide another tool for quantitative analysis. UAS-ST provides 
new capabilities which ASC-U was unable to perform. With the use of the CPLEX 
solver in GAMS, UAS-ST is able to quickly analyze tens of thousands of schedules for 
individual UAVs and then choose a near-optimal subset to provide a complete schedule. 
The user is able to modify all aspects of the UAS structure through a simple spreadsheet 
interface. Output is provided in the same format, which makes understanding the results 
and error tracking much simpler. ASC-U requires the use of several heuristics to meet 
the needs of OUAS. UAS-ST directly addresses these issues without the need for 
additional heuristics. UAS-ST provides a rapid analysis tool and should be implemented 


immediately in the ongoing OUAS study. 
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C. FUTURE DEVELOPMENT 


UAS development is expanding rapidly with several new technologies becoming 
available. Therefore, modeling requirements are expected to change rapidly as well. A 
possible issue with UAS-ST is the requirement for the GAMS CPLEX solver to produce 
a quick solution. Due to portability issues, the future development of a non-proprietary 
heuristic is a valuable topic for future analysis. To improve on UAS-ST, the heuristic 
needs to modify the enumeration routine used in the model. Rather than a total 
enumeration of all possible schedules, a new heuristic should initially select higher value 
routes and not consider low value routes beneath a designated limit. One possible 
method is a greedy heuristic which selects the highest value subset of individual 
schedules and eliminate all others from consideration. By initially eliminating low value 
routes, the overall processing time is greatly reduced. In the long term, analysts in theater 
may run UAS-ST on their laptops as they plan for actual UAS employment. The ultimate 
goal is to find a near-optimal solution without the need for specialty software 


applications. 
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APPENDIX A TEST SCENARIO 


Tier I Mission Areas (m1-m40) 
Tier II Mission Areas (m41-m103) 
Tier III Mission Areas (m104-150) 


Location Mission 


Class 


Value 
Per Time 
Period 


Time 
Periods 


N323934 


W1 144850 





N322947 


W1144606 





N322825 


W1144547 





N322648 
N322549 


W1144730 
W1144914 





N324012 


W1 142426 





N324021 


W1142011 





N324055 


W1142318 





N323922 
N323848 


W1142442 
W1142056 





N325910 


W1142608 





N325733 


W1142738 





N325757 


W1142944 





N325927 
N330047 


W1 142938 
W1142630 





N323934 


W1 144850 





N322947 


W1144606 





N322825 


W1144547 





N322648 
N322549 


W1 144730 
W1144914 





N324012 


W1 142426 





N324021 


W1142011 





N324055 
N323922 


W1142318 
W1142442 





N323848 


W1142056 





N325910 


W1 142608 





N325733 


W1142738 





N325757 


W1142944 





N325927 


W1142938 





N330047 


W1 142630 





N323934 


W1144850 





m33 


N324012 
N325910 


W1 142426 
W1142608 


AV_CLR | 


96 








m34 





N322947 





W1144606 








ISR 





70 











ad 









































































































































m35 N324021 |W1142011 ISR 69 76 8 10 
m36 N325733 |W1142738 ISR 71 82 12 10 
m37 N322825 |W1144547 ISR 73 76 4 10 
m38s N324055 |W1142318 ISR 75 82 8 10 
m39 N325757 |W1142944 ISR 77 88 12 10 
m40 N322648 |W1144730 ISR 79 82 4 10 
m41 N321813 |W1145516 ISR 1 4 4 iS) 
m42 N321817 |W1143321 ISR 3 18 16 15 
m43 N322752 |W1150128 ISR 5 36 32 15 
m44 N322413 |W1140639| AV_CLR 1 48 48 10 
m45 N324657 |W1135856 CR 1 48 48 10 
m46 N325149 |W1142341| LSR_P | 1 8 8 30 
m47 N330310 |W1144304 ISR 7 10 4 15 
m48 N331121 |W1141606 ISR S) 24 16 15 
m49g N331641 |W1144635 ISR 11 42 32 15 
m50 N321813 |W1145516| AV_CLR ) 56 48 10 
m51 N321817 |W1143321 CR S) 56 48 10 
m52 N322752 |W1150128| LSR_P 17 24 8 30 
m53 N322413 |W1140639 ISR 13 16 4 15 
m54 N324657 |W1135856 ISR 15 30 16 15 
m55 N325149 |W1142341 ISR 17 48 32 15 

N330310 |W1144304 

N331121 |W1141606 

N331641 |W1144635 

N321813 |W1145516 

N321817 |W1143321 

N322752 |W1150128 

N322413 |W1140639 

N324657 |W1135856 

N325149 |W1142341 

N330310 |W1144304 
m66 N331121 |W1141606 ISR 27 42 16 1S 
m67 N331641 |W1144635 ISR 29 60 32 15 
m68 N321813 |W1145516| AV_CLR 33 80 48 10 
m69 N321817 |W1143321 CR 33 80 48 10 
m70 N322752 |W1150128| LSR_P 57 64 8 30 
m71 N322413 |W1140639 ISR 31 34 4 LS) 
m72 N324657 |W1135856 ISR 33 48 16 15 
m73 N325149 |W1142341 ISR 35 66 32 15 
m74 N330310 |W1144304) AV_CLR 41 88 48 10 
m75 N331121 |W1141606 CR 41 88 48 10 
m76 N331641 |W1144635| LSR_P 73 80 8 30 
m77 N321813 |W1145516 ISR 37 40 4 15 
m78 N321817 |W1143321 ISR 39 54 16 15 
m79 N322752 |W1150128 ISR 41 72 32 15 
ms0 N322413 |W1140639| AV_CLR 49 96 48 10 
msi N324657 |W1135856 CR 49 96 48 10 
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m8&2 


N325149 


W1142341 


LSR_P 





m83 


N330310 


W1144304 


ISR 





m84 


N331121 


W1141606 


ISR 





m85 


N331641 


W1144635 


ISR 





m86 


N321813 


W1145516 


AV_CLR 





m87 


N321817 


W1143321 


CR 





m88& 


N322752 
N322413 


W1150128 
W1140639 


LSR_P_ 





N324657 


W1135856 





N325149 


W1142341 





N330310 


W1144304 





N331121 
N331641 


W1141606 
W1144635 





N321813 


W1145516 





N321817 


W1143321 





N322752 


W1150128 





N322413 
N324657 


W1140639 
W1135856 





N325149 


W1142341 





N330310 


W1 144304 





N331121 


W1141606 





N331641 
N312538 


W1144635 
W1145510 





N334259 


W1155906 





N333540 


W1115204 





N312538 


W1145510 





N334259 


N333540 


W1155906 
W1115204 





N312538 


W1145510 





N334259 


W1155906 





m113 N312538 |W1145510 CR 96 96 5 


N333540 


N334259 


W1115204 


W1 155906 


SIGINT 


a es 





N333540 


W1115204 


SIGINT 





N312538 


W1145510 


SIGINT 





N334259 
N333540 


W1155906 
W1115204 


SIGINT 





N312538 


W1145510 


o;a/Aa|/a}/a/oa 





N334259 


W1155906 





N333540 


W1115204 





N312538 


W1145510 





N334259 


W1155906 





N333540 


W1115204 





N312538 


W1145510 





N334259 


W1155906 





mi28 





N333540 
N312538 





W1115204 
W1145510 





STK 






































mi29 N334259 |W1155906; STK 49 64 16 60 
m130 N333540 |W1115204| STK 65 80 16 60 
m131 N312538 |W1145510| STK 81 96 16 60 
m132 N334259 |W1155906| STK 25 40 16 60 
m133 N333540 |W1115204| STK 73 88 16 60 
m134 N312538 |W1145510| LSR_D 1 8 8 40 
m135 N334259 [W1155906| LSR_D | 25 32 8 40 


N333540 


W1115204 





N312538 


W1145510 





N334259 


W1155906 





N333540 


W1115204 





N312538 
N334259 


W1145510 
W1155906 





N333540 


W1115204 





N312538 


W1145510 





N334259 


W1155906 





N333540 
N312538 


W1115204 
W1145510 





N334259 


W1155906 





N333540 


W1115204 





N312538 


W1145510 








N334259 





W1155906 
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APPENDIX B VBA CODE 


Sub SolveUAS () 


Dim bSolveResult As Boolean 








bSolveResult = SolveProblem (True) 


End Sub 





Sub BuildData() 
Dim uasDataFN As Integer 
Dim uasSetsFN As Integer 
Dim uasDynSetsFN As Integer 
Dim logFN As Integer 
Dim rTiers As Range 
Dim rGCSs As Range 
Dim rLRSs As Range 
Dim rUAVs As Range 
Dim rMissions As Range 
Dim rPayloads As Range 
Dim rSchedules As Range 
Dim rClasses As Range 


Dim lNumTiers As Long, 1lNumGCSs As Long, l1NumLRSs As Long, 1NumUAVs As 
Long 
Dim 1NumMissions As Long, lNumPayloads As Long, lNumClasses As Long 


Dim 1NumSchedules As Long 








1 

1 
Dim 1NumLRSSchedules As Long 
Dim bNextRoute As Boolean 
Dim rValues As Range 
Dim 1NumValues As Long 
Dim rDependencies As Range 
Dim 1PlanPeriods As Long 
Dim 1FirstPeriod As Long 
Dim lLastPeriod As Long 


Dim 1PeriodsPerHour As Long 





Dim 1Count As Long 
Dim currCell As Range 
Dim iRow As Integer, iCol As Integer 


Dim 1Mission As Long, l1Mission2 As Long, l1UAV As Long 
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' 
assigned} 


Dim 1Period As Long, 1Periodl As Long, 1Period2 As Long 


Dim 1Tier As Long, lPayload As Long, 1GCS As Long, 


Dim 1Class As Long 


Dim tierNames As StringDictionaryClass 
Dim uavNames As StringDictionaryClass 

Dim gcsNames As StringDictionaryClass 

Dim lrsNames As StringDictionaryClass 

Dim missionNames As StringDictionaryClass 
Dim payloadNames As StringDictionaryClass 
Dim classNames As StringDictionaryClass 
Dim periodNames As StringDictionaryClass 


Dim tValues As TupleDictionaryClass 


uasDataFN = FreeFile() 


1LRS As Long 





Open ThisWorkbook.path & "\uasData.gms" For Output As uasDataFN 





uasSetsFN = FreeFile() 


Open ThisWorkbook.path & "\uasSets.gms" For Output As uasSetsFN 


logFN = FreeFile() 


Open ThisWorkbook.path & "\uasXL.log" For Output As logFN 














{STATIC, GROUNDED} SETS 











Vv Tier levels 

h mission types {active, 
g Ground control stations 

1 Launch and recovery sites 

m Missions 

t Time periods 

Pp Payload module types 

Ss UAV employment schedules 
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passive} 


{missions that ships can be 





Print #logFN, "Starting output" 

Print #uasSetsFN, " OPTIONS" 

Print #uasSetsFN, " MIP = " & wsDashboard.Range ("SOLVER") 
Print #uasSetsFN, " OPTCR = " & wsDashboard.Range ("OPTCR") 
Print #uasSetsFN, " RESLIM = " & wsDashboard.Range ("RESLIM") 
Print #uasSetsFN, " ;" 





Print #uasSetsFN, " SETS" 








Set rTiers = wsTiers.Range("A4", wsTiers.Range("A4") .End(x1lDown) ) 





Call GAMSWriteSetDef(rTiers, "v", uasSetsFN, "Tier levels") 


Set rUAVs = wsUAV.Range("A4", wsUAV.Range("A4") .End(x1lDown) ) 





Print #uasSetsFN, " Ay fe 

Print #uasSetsFN, " h_active" 
Print #uasSetsFN, " h_passive" 
Print #uasSetsFN, " i 


Set rGCSs = wsGCS.Range("A4", wsGCS.Range("A4") .End(x1lDown) ) 





Call GAMSWriteSetDef(rGCSs, "g", uasSetsFN, "Ground control stations") 





Set rLRSs = wsLRS.Range("A4", wsLRS.Range("A4") .End(x1lDown) ) 


Call GAMSWriteSetDef(rLRSs, me uasSetsFN, "Launch and = recovery 
sites") 

Set rMissions = wsMissions.Range("A4", 
wsMissions.Range("A4") .End(x1lDown) ) 





Call GAMSWriteSetDef(rMissions, "m", uasSetsFN, "Missions") 


Set rPayloads = wsPayloads.Range("A4", 
wsPayloads.Range ("A4") .End(x1lDown) ) 








' Call GAMSWriteSetDef(rPayloads, "p", uasSetsFN, "Payloads") 


Set rClasses = wsMClasses.Range("A4", 
wsMClasses.Range ("A4") .End(x1lDown) ) 





Call GAMSWriteSetDef(rClasses, "c", uasSetsFN, "Payload Classes") 


1PeriodsPerHour = CLng(wsDashboard.Range ("PERIODS_PER_HOUR") ) 








1PlanPeriods = CLng(wsDashboard.Range ("HORIZON") ) 





lFirstPeriod = CLng(wsDashboard.Range("FIRST_PERIOD") ) 


lLastPeriod = CLng(wsDashboard.Range ("LAST_PERIOD") ) 





If lLastPeriod > 1PlanPeriods Then 


MsgBox "Last Period comes after maximum planning horizon: 
truncating." 
lLastPeriod = 1PlanPeriods 
End If 





If 1FirstPeriod > lLastPeriod Then 





MsgBox "Empty planning horizon: first Period occurs after last 
Period. Aborting." 





GoTo subExit 





End If 

Print #logFN, "Periods " & 1PlanPeriods & " First: " & 1lFirstPeriod & 
"Last: " & lLastPeriod 

Print #uasSetsFN, " t Time Periods /t_" & Format(lFirstPeriod) & 


"xt" & Format (lLastPeriod) & "/" 
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Set tierNames = New StringDictionaryClass 





If DefineGroundedSet (rTiers, tierNames) = False Then 





Print #logFN, "Duplicate element in grounded set (tiers). User chose 
to cancel: no further output." 








GoTo subExit 


End If 





lNumTiers = tierNames.Count 








For 1Count = 1 To 1NumTiers 
Print #logFN, Format (1lCount) & " : " & tierNames.Item(1Count) 


Next 1Count 


Set uavNames = New StringDictionaryClass 


If DefineGroundedSet (rUAVs, uavNames) = False Then 








Print #logFN, "Duplicate element in grounded set (uavs). User chose 
to cancel: no further output." 





GoTo subExit 


End If 





1NumUAVs = uavNames.Count 


For 1Count = 1 To 1NumUAVs 
Print #logFN, Format (lCount) & " : " & uavNames.Item(1Count) 


Next 1Count 


Set missionNames = New StringDictionaryClass 


If DefineGroundedSet (rMissions, missionNames) = False Then 








Print #logFN, "Duplicate element in grounded set (missions). User 
chose to cancel: no further output." 





GoTo subExit 


End If 





1NumMissions = missionNames.Count 
For 1Count = 1 To 1NumMissions 
Print #logFN, Format (1lCount) & ": " & missionNames.Item(1Count) 


Next 1Count 


Set payloadNames = New StringDictionaryClass 








If DefineGroundedSet (rPayloads, payloadNames) = False Then 


Print #logFN, "Duplicate element in grounded set (payloads). User 
chose to cancel: no further output." 
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GoTo subExit 


End If 





lNumPayloads = payloadNames.Count 


For 1lCount = 1 To 1lNumPayloads 
Print #logFN, Format (1lCount) & " : " & payloadNames.Item(1Count) 


Next 1Count 


Set classNames = New StringDictionaryClass 


If DefineGroundedSet (rClasses, classNames) = False Then 








Print #logFN, "Duplicate element in grounded set (classes). User 
chose to cancel: no further output." 





GoTo subExit 


End If 





lNumClasses = classNames.Count 


For 1lCount = 1 To 1NumPayloads 
Print #logFN, Format(1lCount) & " : " & payloadNames.Item(1Count) 


Next 1Count 


Set gcsNames = New StringDictionaryClass 


If DefineGroundedSet (rGCSs, gcsNames) = False Then 








Print #logFN, "Duplicate element in grounded set (GCS). User chose to 
cancel: no further output." 





GoTo subExit 


End If 





1NumGCSs = gcsNames.Count 


For 1Count = 1 To 1NumGCSs 
Print #logFN, Format (l1Count) & " : " & gcsNames.Item(1Count) 


Next 1Count 


Set lrsNames = New StringDictionaryClass 


If DefineGroundedSet (rLRSs, lrsNames) = False Then 








Print #logFN, "Duplicate element in grounded set (LRS). User chose to 
cancel: no further output." 





GoTo subExit 


End If 





1NumLRSs = lrsNames.Count 
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For 1Count = 1 To 1NumLRSs 
Print #logFN, Format (lCount) & " : " & IlrsNames.Item(1Count) 


Next 1Count 


Set periodNames = New StringDictionaryClass 

periodNames.Size = 3 * 1PlanPeriods 

For 1lCount = 1 To 1PlanPeriods 
periodNames.Add "t" & Format (1Count) 


Next 1Count 


For 1Count = 1FirstPeriod To lLastPeriod 
Print #logFN, Format(lCount) & " : " & periodNames.Item(1Count) 


Next 1Count 





' {GIVEN} PARAMETERS 


























i num_uavs(l,v,t) number of uavs available at lrs 1 in tier v in period 
G 
num_mods (p,t) number of payload modules available of ID p in period 
& 
‘ gcs_cap(g,t,h) max num UAVs on gcs g in period t on mission type h 
lrs_cap(1,t) max num UAVs at LRS 1 in period t 
: length (m) maximum periods of coverage of mission m 
‘ val (m) value of mission m 
Y k(1,s,t) indicator if lrs 1 used by schedule s in period t 
a 
; 
Print #uasDataFN, " PARAMETERS" 
i num_uavs (1,v,t) number of uavs available at lrs 1 in tier v in 
period t 
Print #uasDataFN, " num_uavs(l,v,t) /" 
ReDim l1NumUAVinPeriod(1l To lrsNames.Count, aL To tierNames.Count, 


lFirstPeriod To lLastPeriod) As Long 
For 1Period = 1FirstPeriod To lLastPeriod 
For 1LRS = 1 To lrsNames.Count 
For lTier = 1 To tierNames.Count 


1NumUAVinPeriod(1LRS, lTier, 1Period) = 0 





Next lTier 
Next 1LRS 
For iRow = 1 To rUAVs.Rows.Count 


If rUAVs(iRow, COL_UAV_FIRST_AVAIL) <= 1Period And rUAVs(iRow, 
COL_UAV_LAST_AVAIL) >= 1Period Then 
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1LRS = lrsNames.Lookup (rUAVs (iRow, COL_UAV_LRS) ) 


1lTier = tierNames.Lookup (rUAVs (iRow, COL_UAV_TIER) ) 





If 1LRS = 0 Then 
MsgBox "Error: invalid LRS name in UAV list" 
wsUAV.Activate 
rUAVs (iRow, COL_UAV_LRS) .Select 


GoTo subExit 





ElseIf 1Tier = 0 Then 





MsgBox "Error: invalid tier name in UAV list" 


wsUAV.Activate 


rUAVs (iRow, COL_UAV_TIER) .Select 








GoTo subExit 





Else 


1NumUAVinPeriod(1LRS, 1Tier, 1Period) = I1NumUAVinPeriod(1LRS, 
lTier, 1Period) + 1 


End If 





End If 





Next iRow 


Next 1Period 


For 1LRS = 1 To lrsNames.Count 
For lTier = 1 To 1NumTiers 
For 1Period = 1FirstPeriod To lLastPeriod 
If 1lNumUAVinPeriod(1LRS, l1Tier, 1Period) >= 1 Then 


Print #uasDataFN, " 1." & lrsNames.Item(1LRS) & ".v_" & 
tierNames.Item(lTier) & ".t_" & 1Period & "" & 





1NumUAVinPeriod(1LRS, l1Tier, 1Period) 
End If 





Next 1Period 
Next lTier 
Next 1LRS 


Print #uasDataFN, " rhs 


ReDim numMods(classNames.Count, tierNames.Count) As Long 
For 1Payload = 1 To rPayloads.Rows.Count 
iRow = classNames.Lookup(rPayloads(lPayload, COL_PAYLOAD_CLASS) ) 


iCol = tierNames.Lookup(rPayloads(lPayload, COL_PAYLOAD_TIER) ) 





numMods (iRow, iCol) = numMods (iRow, iCol) + rPayloads (1Payload, 
COL_PAYLOAD_QUANTITY) 


Next 1lPayload 
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num_mods (c,v) number of payload modules available of class c for 


tier v 
Print #uasDataFN, " num_mods(c,v) /" 
For 1Class = 1 To classNames.Count 
For lTier = 1 To tierNames.Count 
Print #uasDataFN, uy aes & classNames.Item(1Class) & ‘ies | eee & 
tierNames.Item(lTier) & " " & numMods(l1Class, l1Tier) 





Next lTier 
Next 1Class 


Print #uasDataFN, " ae 


gcs_cap(g,t,h) max num UAVs on gcs g in period t on mission type h 
Print #uasDataFN, " gcs_cap(g,t,h) /" 
For 1GCS = 1 To 1NumGCSs 

For 1Period = 1FirstPeriod To lLastPeriod 


Print #uasDataFN, " g_" & gcsNames.Item(1GCS) & ".t_" & 1Period 
& ".h_active " & rGCSs(1GCS, COL_GCS_ACTIVE) 





Print #uasDataFN, " g_" & gcsNames.Item(1GCS) & ".t_" & 1lPeriod 
& ".h_passive " & rGCSs(1GCS, COL_GCS_PASSIVE) 





Next 1Period 


Next 1GCS 
Print #uasDataFN, " au 

‘ lrs_cap(1,t) max num UAVs at LRS 1 in period t 
Print #uasDataFN, " lrs_cap(1,t) /" 


For 1LRS = 1 To 1NumLRSs 
For 1Period = 1FirstPeriod To lLastPeriod 


Print #uasDataFN, " 1_" & lrsNames.Item(1LRS) & ".t_" & 1Period 
& " " & rLRSs(1LRS, COL_LRS_ACTIVE) 








Next 1Period 








Next LLRS 
Print #uasDataFN, " as 
: length (m) maximum periods of coverage of mission m 
Print #uasDataFN, " length(m) /" 
For 1Mission = 1 To 1NumMissions 
Print #uasDataFN, " m_" & missionNames.Item(1lMission) & " " & 

rMissions(1lMission, COL _MISSION_PERIODS) 

Next 1Mission 


Print #uasDataFN, " is 
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: val (m) value of mission m 


Print #uasDataFN, " val(m) /" 
For 1Mission = 1 To l1NumMissions 
Print #uasDataFN, " m_" & missionNames.Item(1lMission) & " " & 


rMissions(lMission, COL_MISSION_VALUE) 





Next 1Mission 


Print #uasDataFN, " ia 





Calculate distances between GCS-Mission pairs, LRS-mission pairs, and 
mission-mission pairs 


ReDim dDistGM(lNumGCSs, lNumMissions) As Double 
ReDim dDistLM(1lNumLRSs, 1NumMissions) As Double 


ReDim dDistMM(1lNumMissions, 1lNumMissions) As Double 


Dim i As Long, j As Long 


For 1GCS = 1 To rGCSs.Rows.Count 


i = gcsNames.Lookup(rGCSs(1GCS, COL_GCS_NAME) ) 





For 1Mission = 1 To rMissions.Rows.Count 





) 


dDistGM(i, 5) = SphericalDistance(rGCSs (1GCS, COL_GCS_LAT), 
rGCSs(1GCS, COL_GCS_LON), 


j = missionNames.Lookup(rMissions(lMission, COL_MISSION_NAME 


rMissions (1Mission, 
COL_MISSION_LAT), rMissions(l1Mission, COL_MISSION_LON) ) 


Next 1Mission 


Next 1GCS 


For 1LRS = 1 To rLRSs.Rows.Count 





i = lrsNames.Lookup(rLRSs(1LRS, COL_LRS_NAME) ) 


For 1Mission = 1 To rMissions.Rows.Count 





) 


dDistLM(i, 5) = SphericalDistance(rLRSs (1LRS, COL_LRS_LAT), 
rLRSs(1LRS, COL_LRS_LON), 


j = missionNames.Lookup(rMissions(1lMission, COL_MISSION_NAME 





rMissions (1Mission, 
COL_MISSION_LAT), rMissions(l1Mission, COL_MISSION_LON) ) 


Next 1Mission 


Next 1LRS 


For 1Mission = 1 To rMissions.Rows.Count 





i = missionNames. Lookup (rMissions(1Mission, COL_MISSION_NAME) ) 


For 1Mission2 = 1 To rMissions.Rows.Count 
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j = missionNames.Lookup(rMissions(lMission2, COL_MISSION_NAME) ) 





dDistMM (i, 5) = SphericalDistance(rMissions (1Mission, 
COL_MISSION_LAT), rMissions(1lMission, COL_MISSION_LON), 


rMissions(l1Mission2, 
COL_MISSION_LAT), rMissions(l1Mission2, COL _MISSION_LON) ) 





Next 1Mission2 


Next 1Mission 


"Find GCSs within range of each mission. If none, requir on of the 
AV_CLR missions 





Dim pFN As Integer, hpFN As Integer 
Dim numPreq As Long 

numPreq = 0 

pFN = FreeFile() 

Open "p.gms" For Output As pFN 

hpFN = FreeFile() 

Open "hp.gms" For Output As hpFN 


Print #pFN, " SET P(m,mp) /" 





Print #hpFN, " SET HP(m) /" 





For 1Mission = 1 To missionNames.Count 
For 1GCS = 1 To gcsNames.Count 


If dDistGM(1GCS, l1Mission) <= rGCSs(1GCS, COL_GCS_LOS) Then GoTo 
nextMissionGCs 


Next 1GCS 
numPreq = numPreq + 1 


Print #hpFN, " m_" & missionNames.Item(1Mission) 





For 1Mission2 = 1 To rMissions.Rows.Count 


Ae rMissions(l1Mission2, COL_MISSION_CLASS) = 
wsGCS.Range ("LINK") .Value Then 


Print #pFN, " m_" & missionNames.Item(1lMission) & ".m_" & 


rMissions(1lMission2, COL_MISSION_NAME 





End If 





Next 1Mission2 
nextMissionGCs: 
Next 1Mission 
If numPreq = O Then 


Print #pFN, mn m_" &  missionNames.Item(1) & mS" & 
missionNames.Item(1) 


Print #hpFN, " m_" & missionNames.Item(1) 
End If 


Print #pFN, " ‘aa 








Print #pFN, " ;" 
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gcsNames. 


Close pFN 
Print #hpFN, " fe 
Print #hpFN, " ;" 


Close hpFN 


Dim mgFN As Integer 


mgFN = FreeFile() 


Open "mg.gms" For Output As mgFN 





Print #mgFN, " SET mg(m,g) /" 
For 1Mission = 1 To missionNames.Count 
For iRow = 1 To rGCSs.Rows.Count 


1GCS = gcsNames. Lookup (rGCSs (iRow, 


If dDistGM(1GCS, l1Mission) 


Print #mgFN, " 
Item(1GCS) 


End If 





Next iRow 
Next 1Mission 
Print #mgFN, " /;" 


Close mgFN 


Dim mhFN As Integer 


mhFN = FreeFile() 


m" 


COL_GCS_NAME) ) 





<= rGCSs(iRow, COL_GCS_LOS) Then 


& missionNames.Item(1lMission) & ". 


Open "mh.gms" For Output As mgFN 


ReDim bClassActive(classNames.Count) 


For iRow = 1 To rClasses.Rows.Count 


As Boolean 


1Class = classNames.Lookup(rClasses (iRow, 


If CLng(rClasses (iRow, 








bClassActive(lClass) = True 
Else 
bClassActive(l1Class) = False 
End If 
Next iRow 
Print #mhFN, " SET mh(m,h) /" 
For iRow = 1 To rMissions.Rows.Count 
1Class = classNames.Lookup(rMissions (iRow, 
If bClassActive(l1Class) Then 


COL_MCLASS_ACTIVE) ) 
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COL_MCLASS_NAM 


Gl 


= 1 Then 


COL_MISSION_CLASS) ) 





Print #mhFN, i m_" & rMissions(iRow, COL_MISSION_NAME) & 
" h_active" 





Else 


Print #mhFN, m_" & rMissions(iRow, COL_MISSION_NAME) & 
",h_passive" 





End If 





Next iRow 
Print #mhFN, " /;" 


Close mhFN 


Dim gvFN As Integer 
gvEN = FreeFile() 
Open "gv.gms" For Output As mgFN 


Print #gvFN, " SET GV(g,v) /" 





For iRow = 1 To rGCSs.Rows.Count 


1GCS = gcsNames.Lookup(rGCSs(iRow, COL_GCS_NAMBE) ) 





lTier = tierNames.Lookup(rGCSs (iRow, COL_GCS_TIER) ) 





Print #QVEN, " gas & gcsNames.Item(1GCS) & Wot & 
tierNames.Item(l1Tier) 





Next iRow 
Print #gvFN, " /;" 


Close gvFN 


Print #uasDataFN, " r(m) /" 




















For iRow = 1 To rMissions.Rows.Count 
1Mission = missionNames.Lookup (rMissions (iRow, COL_MISSTION_NAME) ) 
If Trim(CStr(rMissions(iRow, COL_MISSION_REQUIRED))) <> "" Then 
TE CLng (rMissions (iRow, COL_MISSION_REQUIRED) ) > 











CLng (rMissions(iRow, COL_MISSION_PERIODS)) Then 


If MsgBox("Required coverage exceeds total mission time", 
vbOKCancel) = vbCancel Then 


wsMissions.Activate 





rMissions(iRow, COL_MISSION_REQUIRED) .Select 

















Exit Sub 
Else 
Print #uasDataFN, " m_" & missionNames.Item(lMission) & " m 


& CLng(rMissions(iRow, COL_MISSION_PERIODS) ) 





End If 


GI 
Bb 





se 





Print #uasDataFN, " m_" & missionNames.Item(lMission) & " ie 
CLng (rMissions(iRow, COL_MISSION_REQUIRED) ) 














End If 
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' 
period t 
' 
period t 


mission type h 


Else 





End If 


Next iRow 


Print #uasDataFN, " m_" & missionNames.Item(1lMission) & " 0" 


Print #uasDataFN, " /" 


Count number of UAVs in each tier at each LRS 





ReDim tierAtLRS (tierNames.Count, lrsNames.Count) As Long 


For 1LRS = 1 To lrsNames.Count 


For lTier 


= 1 To tierNames.Count 


tierAtLRS(lTier, 1LRS) = 0 


Next lTier 


Next 1LRS 


For 1LUAV = 1 


To rUuAVs.Rows.Count 


i = tierNames. Lookup (rUAVs (LUAV, COL_UAV_TIER) ) 





j = lrsNames.Lookup (rUAVs (1UAV, COL_UAV_LRS) ) 


tierAtLRS(i, j) = tierAtLRS(i, j) + 1 


Next LUAV 


a(m,s,t) 


b(v,s,t) 


f(g,s,t,h) 


k(1,s,t) 


q(p,s,t) 


indicator if mission m covered by schedule s in 


indicator if UAV in Tier v required by schedule s in 


indicator if gcs g used by schedule s in period t for 


indicator if lrs 1 used by schedule s in period t 


indicator if payload p used by schedule s in period t 


Dim 1MaxRoutes As Long 


1MaxRoutes = 








CLng (wsDashboard. Range ("MAX_ROUTES") ) 


Dim 1MaxMissions As Long 


1MaxMissions 


= CLng(wsDashboard.Range ("MAX_MISSIONS")) + 1 


ReDim a(lMaxRoutes, wsDashboard.Range("LAST_PERIOD")) As Long 


ReDim b(1MaxRoutes) As Long 




















ReDim f_act (lMaxRoutes, wsDashboard.Range("LAST_PERIOD")) As Long 
ReDim f_pas(lMaxRoutes, wsDashboard.Range("LAST_PERIOD")) As Long 
ReDim k(1lMaxRoutes, wsDashboard.Range("LAST_PERIOD")) As Long 
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COL_MISSION_NAME 


COL_MISSION_NAME 


ReDim q(lMaxRoutes, classNames.Count) 





ReDim stack(lMaxMissions + 1) As Long 
ReDim onStack(lNumMissions) As Long 
ReDim dwell(lMaxMissions + 1) As Long 


For 1Mission 1 To 1NumMissions 














As 


Long 


onStack(lMission) = 0 
Next 1Mission 
Dim top As Long 
top = 0 
lNumSchedules = 0 
1NumLRSSchedules = 0 
bNextRoute = Tru 
1LRS = 1 
Dim s As Long, t As Long 
Dim 1NextMission As Long, l1CurrMission As Long 
ReDim 1NextStack(lMaxMissions + 1) As Long 
Dim dRemTime As Double 
Dim 1RemPeriods As Long 
Dim dTransitTime As Double 
Dim lTransitPeriods As Long 
Dim lcurrNumClasses As Long 
Dim dTierSpeed As Double 
Dim 1ThisLaunchPeriod As Long 
Dim 1ThisRecoverPeriod As Long 
Dim 1MaxClasses As Long 
ReDim lTransitPeriod(lMaxMissions + 1) As Long 
ReDim lIdlePeriod(1MaxMissions + 1) As Long 
ReDim 1MissionPeriod(lMaxMissions + 1) As Long 
ReDim 1StartPeriod(lNumMissions) As Long 
ReDim 1FinishPeriod(lNumMissions) As Long 


For 1Mission 


1 To rMissions.Rows.Count 


1StartPeriod(missionNames. Lookup (rMissions (1Mission, 


CLng (rMissions (1lMission, 


Col 


L MISSION _START_P 


1lFinishPeriod(missionNames. Lookup (rMissions(1lMission, 








CLng (rMissions (1lMission, 


Next 1Mission 
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Col 





L MISSION_FINISH FE 


Dp 





ERIOD) ) 


ERIOD) ) 





ReDim dMissionValue (missionNames.Count) 


For 1Mission = 1 To rMissions.Rows.Count 


As Double 


dMissionValue (missionNames.Lookup(rMissions (1Mission, 


COL_MISSION_NAME))) = CDbl1(rMissions(1Mission, 


Next 1Mission 


ReDim lHasClass(lNumClasses) As Long 





For 1Class = 1 To 1NumClasses 
lHasClass(lClass) = 0 

Next 1Class 

lcurrNumClasses = 0 


ReDim reqClass (1NumMissions) 


COL_MISSTION_VALUE 


For 1Mission = 1 To rMissions.Rows.Count 


reqClass (missionNames.Lookup (rMissions (1Mission, 


classNames.Lookup(rMissions (1Mission, 


Next 1Mission 


ReDim 1ClassTierAvail (classNames.Count, 


For 1Class = 1 To classNames.Count 

For lTier = 1 To tierNames.Count 

1ClassTierAvail(l1Class, l1Tier) 
Next lTier 


Next 1Class 


For 1Payload = 1 To rPayloads.Rows.Count 





) 





COL_MISSION_NAM 


COL_MISSION_CLASS) ) 


tierNames.Count) As Long 


1Class = classNames.Lookup (rPayloads (1lPayload, COL_PAYLOAD_CLASS) ) 





lTier = tierNames.Lookup(rPayloads(lPayload, COL_PAYLOAD_TIER) ) 


1ClassTierAvail(1Class, lTier) 


COL_PAYLOAD_QUANTITY) ) 


Next 1lPayload 


Dim aFN As Integer, bFN As Integer, 


As Integer 
aFN = FreeFile() 
Open ThisWorkbook.path & "\a.gms" 
bFN = FreeFile() 
Open ThisWorkbook.path & "\b.gms" 
kFN = FreeFile() 


Open ThisWorkbook.path & "\k.gms" 











QFN = FreeFile() 




















Open ThisWorkbook.path & "\q.gms" 


s = 0 


Print #aFN, " SET A(s,m,t) /" 
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For 


For 


For 


For 


Ou 


Ou 


Ou 


Ou 


tput 


tput 


tput 





tput 


As 


As 


As 


As 


LF. | 


CLng (rPayloads (1Payload, 


fFN As Integer, 


aF 


bF 


kF 


qr 








kFN As Integer, 


qrN 





T, 


Print #bFN, " SET B(s,v,t) / 











Print #qFN, " SET Q(s,c,t) / 


Dim tierLRScombos As Long 


tierLRScombos = 0 


For lTier = 1 To tierNames.Count 


For 1LRS = 1 To lrsNames.Count 


If tierAtLRS(lTier, 1LRS 


) Then 


tierLRScombos = tierLRScombos + 1 


End If 





Next 1LLRS 
Next lTier 
Dim currCombo As Long 


currCombo = 1 


For 1LRS = 1 To rLRSs.Rows.Count 








Print #LogFN, "LRS i & rLRSs(1LRS, COL_LRS_NAME) & meen & 
lrsNames.Lookup(rLRSs(1LRS, COL_LRS_NAME) ) 
For lTier = 1 To rTiers.Rows.Count 
Print #logFN, " tier " & rTiers(lTier, COL_TIER_NAME) & ":" & 














tierNames.Lookup(rTiers(lTier, COL_TI 





ER_NAME) ) 








EE tierAtLRS (tierNames. Lookup (rTiers (1Tier, 


lrsNames.Lookup(rLRSs(1LRS, COL_LRS_N 


AME))) > O Then 














dTierSpeed = rTiers(lTier, COL_TIER_SPEED) 




















dRemTim = xrTiers(lTier, COL_TIER_ENDURANC 





COL_TIER_SPEED) 

















lRemPeriods = Ceiling(dRemTime * 1PeriodsPerH 














ia) 








Q 
(e) 
iE 
Hi 





ER_NAME)), 








'x* rTiers(lTier, 





our) 


1MaxClasses = rTiers(lTier, COL_TIER_CAPACITY) 





Print #logFN, " tierAtLRS? yes 1MaxClasses=" 


For 1Class = 1 To 1NumClasses 


lHasClass(lClass) = 0 


Next 1Class 














lcurrNumClasses = 0 
top = 1 
stack(top) = EVENT_LR 'L/R event at 1LRS 
1NextStack (top) = 1 
Do 
1CurrMission = stack (top) 


1NextMission = 1NextStack (top) 


"Find next mission to put on the stack 
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& 1MaxClasses 


Do While 1NextMission <= 1NumMissions 
If onStack(lNextMission) > 0O Then 
"Print #logFN, " onstack" 
GoTo skipMission 


End If 





Lf dDistLM(1LRS, 1NextMission) > rTiers(lTier, 


COL_TIER_RADIUS) Then 





"Print #logFN, "Out of range" 
GoTo skipMission 


End If 





If 1ClassTierAvail (reqClass(lNextMission), lTier) = 0 Then 
"Print #logFN, "No Payloads available for this Tier" 
GoTo skipMission 


End If 








If 1CurrMission = EVENT_LR Then 











dTransitTime = dDistLM(1LRS, 1NextMission) / dTierSpeed 











Else 
dTransitTime = dDistMM(1CurrMission, 1NextMission) 
dTierSpeed 
End If 
lTransitPeriods = Ceiling(dTransitTime * 1PeriodsPerHour) 





If 1CurrMission <= 0 Then 


nm 1FinishPeriod(1lNextMission) 
1StartPeriod(lCurrMission) + l1TransitPeriods Then 


"Print #logFN, " time warp" 
GoTo skipMission 


End If 








Elself dMissionValue (1NextMission) 
dMissionValue(lCurrMission) Then 


Tf 1FinishPeriod(lNextMission) 
1FinishPeriod(1lCurrMission) + lTransitPeriods Then 


"Print #logFN, " time warp" 
GoTo skipMission 


End If 








Else 


If 1StartPeriod(lNextMission) <= 1StartPeriod(lCurrMission) 


+ lTransitPeriods Then 
"Print #logFN, "  time-value warp" 
GoTo skipMission 


End If 





End If 
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If lTransitPeriods + Ceiling(lPeriodsPerHour * (dDistLM(1LRS, 
1NextMission) / rTiers(lTier, COL_TIER_SPEED))) > 1lRemPeriods Then 

















"Print #logFN, " out of time" 


GoTo skipMission 





End If 
Te (lcurrNumClasses = 1MaxClasses And 
lHasClass(reqClass(lNextMission)) = 0) Then 
"Print #logFN, " no payload space left" 


GoTo skipMission 


End If 





GoTo foundMission 
skipMission: 
1NextMission = 1NextMission + 1 
Loop 


foundMission: 

















If 1NextMission > 1NumMissions Or top >= 1MaxMissions Then ‘out 
of missions: pop stack 
If 1lCurrMission > 0 Then 
lHasClass (reqClass (1CurrMission) ) = 
lHasClass(reqClass(lCurrMission)) - 1 
If lHasClass(reqClass(lCurrMission)) = 0 Then 
lcurrNumClasses = lcurrNumClasses - 1 
End If 
onStack(lCurrMission) = 0 
End If 
top = top - 1 
GoTo nextMission 
End If 
If 1lCurrMission = EVENT_LR Then 
lTransitPeriod(top) = 1TransitPeriods 
1ThisLaunchPeriod = 1StartPeriod(lNextMission) 
lTransitPeriods - 1 
If 1ThisLaunchPeriod < 1 Then 
1ThisLaunchPeriod = 1 
End If 
1RemPeriods = Ceiling (dRemTime x 1PeriodsPerHour) + 
l1ThisLaunchPeriod - 1 
lTransitPeriod(top) = 1ThisLaunchPeriod + l1TransitPeriods 
lIidlePeriod(top) = 1ThisLaunchPeriod - 1 
Else 
1MissionPeriod(top) = 1FinishPeriod(lCurrMission) 
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If dMissionValue(lCurrMission) < dMissionValue(1NextMission) 
And 


1StartPeriod(l1NextMission) - lTransitPeriods - 1 < 
1FinishPeriod(lCurrMission) Then 





1MissionPeriod (top) = lStartPeriod(1NextMission) - 
lTransitPeriods - 1 
End If 
lTransitPeriod(top) = 1MissionPeriod(top) + 1lTransitPeriods 
lIidlePeriod(top) = 1StartPeriod(1NextMission) - 1 
End If 
1NextStack (top) = 1NextMission + 1 
top = top + 1 




















stack(top) = 1NextMission 

onStack(1lNextMission) = top 

1MissionPeriod(top) = 1FinishPeriod(lNextMission) 

1TransitPeriod (top) = 1FinishPeriod (top) + 
Ceiling ((dDistLM(1LRS, 1NextMission) vi rTiers(lTier, COL_TIER_SPEED) ) * 
1PeriodsPerHour) 

1NextStack (top) = 1 

If lHasClass(reqClass(1NextMission)) = 0 Then 

lcurrNumClasses = lcurrNumClasses + 1 
End If 





lHasClass (reqClass (1NextMission) ) = 
lHasClass (reqClass(lNextMission)) + 1 


"New feasible schedule: now copy stack info into appropriate 





arrays 
s=sH+1 
t= l 


"calculate start period of mission sequence: find start time of 
first mission 


" in stack(2), figure out latest departure period to arrive at 
first mission (if 





' earlier than period 1, then 1/r in period 1, transit, and do 
as much of mission 


"as possible 


Print #logEN, "dTransitTime " & dTransitTime & " 
lTransitPeriods " & lTransitPeriods 
Print #logFN, "s:" & CStr(s) & " "; 


For iRow = 1 To top 
Print #logFN, stack(iRow) & ","; 
Next iRow 


Print #logFN, "0" 
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Print #logFN, "m:" & CStr(s) & " "; 
For iRow = 1 To top 
Print #logFN, 1MissionPeriod(iRow) & ","; 
Next iRow 
Print #logFN, "0" 
Print #logFN, "t:" & CStr(s) & " "; 
For iRow = 1 To top 


Print #logFN, lTransitPeriod(iRow) & ","; 





Next iRow 

Print #logFN, "0" 

Print #logFN, "i:" & CStr(s) & " "; 
For iRow = 1 To top 


Print #logFN, lIdlePeriod(iRow) & ","; 





Next iRow 


Print #logFN, "0" 





Do While t < 1ThisLaunchPeriod And t <= lLastPeriod 

















a(s, t) = EVENT_IDLE 





a(s, t) = EVENT_LR 











t=ttil 


Do While t <= 1TransitPeriod(1) And t <= lLastPeriod 





a(s, t) = EVENT_TRANSIT 











f_pas(s, t) =1 
SG 
Loop 


For iRow = 2 To top 
1CurrMission = stack (iRow) 


If iRow < top Then 











1NextMission = stack(iRow + 1) 
Else 

1NextMission = EVENT_LR 'L/R event 
End If 








'" XXX Need upper bounds: time horizon, flight time, etc... 


Do While t <= 1MissionPeriod(iRow) And t <= lLastPeriod 


"fill in arrays 
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f_act(s, t) =1 
t=tet il 
Loop 


Do While t <= 1TransitPeriod(iRow) And t <= lLastPeriod 














a(s, t) = EVENT_TRANSIT 
f_pas(s, t) = 1 
t=tt+il 

Loop 


Do While t <= lIdlePeriod(iRow) And t <= lLastPeriod 




















a(s, t) = EVENT_IDLE 
f_pas(s, t) =1 
t=ttl 

Loop 


Next iRow 


If t <= lLastPeriod Then 








a(s, t) = EVENT_LR 








1ThisRecoverPeriod = t 





f_act(s, t) =1 

k(s, t) = 1LRS 

b(s) = l1Tier 
End If 





For iRow = 1 To lLastPeriod 
Print #logFN, ":" & a(s, iRow); 

Next iRow 

Print #logFN, "" 

For iRow = 2 To top 


For 1Period = 1 To lLastPeriod 


If a(s, 1lPeriod) = stack(iRow) Then 
Print #aFN, ut sn & CStr(s) & wom" & 
missionNames.Item(stack(iRow)) & ".t_" & CStr(lPeriod) 
End If 





Next 1Period 
Next iRow 
For 1Period = 1 To lLastPeriod 


Tt 1Period >= 1ThisLaunchPeriod And 1Period <= 
l1ThisRecoverPeriod Then 


Print #bEN, y suf & CStr(s) & we & 
Trim(rTiers(lTier, COL_TIER_NAME)) & ".t_" & CStr(l1Period) 

















End If 





Next 1Period 











Print #kFN, " k('s_" & CStr(s) & "','1l_" & Trim(rLRSs(1LRS, 
COL_LRS_NAME)) & "', 't_" & CStr(lThisLaunchPeriod) & "')=1;" 
Print #kFN, " k('s_" & CStr(s) & "','1l_" & Trim(rLRSs(1LRS, 
COL_LRS_NAME)) & "', 't_" & CStr(lThisRecoverPeriod) & "')=1;" 
For 1Class = 1 To classNames.Count 
If lHasClass(1lClass) > 0 Then 
q(s, 1Class) = 1 
For 1Period = 1 To lLastPeriod 
If 1Period >= 1ThisLaunchPeriod And 1Period <= 
l1ThisRecoverPeriod Then 
Print #qFN, u Ss" -&. <cStr.(s) & Me & 


classNames.Item(l1Class) & ".t_" & CStr(lPeriod) 


End If 





Next 1Period 


End If 





Next 1Class 


If s >= 1MaxRoutes Then GoTo doneMissions 











If s Mod 1000 = 0 Then wsDashboard.Range("ROUTESCELL") = s 
If s >= (currCombo - 1) * l1MaxRoutes / (tierLRScombos) Then 
GoTo nextTier 
nextMission: 


Loop While top > 0 
End If 





nextTier: 
currCombo = currCombo + 1 


Next lTier 
































Next 1LLRS 
doneMissions: 
wsDashboard. Range ("ROUTESCELL") = s 
Print #uasSetsFN, " s /s_1*s_" & CStr(s) & " /" 
Print #uasSetsFN, " ;" 
Print #uasDataFN, " ;" 
Print #aFN, " /" 
Print #aFN, " ;" 
Print #bFN, " /* 
Print #bFN, " ;" 
Print #qFN, " a 
Print #qFN, " ;" 
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subExit: 


Print 
Close 
Close 
Close 
Close 
Close 
Close 
Close 


End Sub 





#logFN, " 
uasSetsFN 
uasDataFN 
aFN 
bFN 
kFN 
qFN 


LogFN 





Ending output" 
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APPENDIX C GAMS CODE 


$TITLE UAS Planner 08.06.02 


$ontext 

MCCDC/OAD UAS Mission Planner 

Assigns a fleet of UAVs to a list of missions with various 
requirements and time windows over a fixed time horizon of many 
time periods. 

Solved using explicit column generation, one column per UAS schedule. 
Sofftext 

$inlinecom { } 

FILE UASP_log/UASP.log/ ; 

PUT UASP_log ; 

PUT 'UASP - MCCDC/OAD UAS Planner 08.06.02' // ; 
PUTCLOSE UASP_log ; 


FILE UASP_sta/UASP-.sta/ ; 
PUT UASP sta ; 

PUT 'ERROR - Did not solve' // ; 
PUTCLOSE UASP_sta ; 


FILE VALUE_CSV/VALUE.csv/ ; 

PUT VALUE_CSV ; 

PUT''/; { scratch legacy file contents } 
PUTCLOSE VALUE_CSV ; 


FILE MISSIONS_CSV/MISSIONS.csv/ ; 
PUT MISSIONS_CSV ; 

PUT''/; { scratch legacy file contents } 
PUTCLOSE MISSIONS_CSV ; 


65 


FILE SCHEDULES_CSV/SCHEDULES.csv/ ; 
PUT SCHEDULES_CSV ; 

PUT''/; { scratch legacy file contents } 
PUTCLOSE SCHEDULES_CSV 


FILE DUALS_CSV/DUALS.csv/ ; 

PUT DUALS_CSV ; 

PU ds { scratch legacy file contents } 
PUTCLOSE DUALS_CSV 


PUT UASP_log ; 
UASP_log.ap=1; _{ re-open log file, appending to prior contents } 


OPTIONS 
optcr =0.1 
limrow = 0 
limcol = 0 
SEED = 1234 


ITERLIM = 10000000 
LP =CPLEX 
RMIP =CPLEX 
MIP =CPLEX 
S$ONTEXT 
{STATIC, GROUNDED} SETS 


Vv Tier levels 

h mission types {active, passive} 

g Ground control stations 

] Launch and recovery sites 

m Missions {missions that ships can be assigned} 
c Payload classes 
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t Time periods 

S UAV emplyment schedules 
$OFFTEXT 
$INCLUDE uasSets.gms 


alias(m,mp); 
S$ONTEXT 
{GIVEN} PARAMETERS 
num_uavs(l,v,t) number of uavs available in tier v in period t 
num_mods(c,t) number of payload modules available of class c in period t 


gcs_cap(g,t,h) | max num UAVs on gcs g in period t on mission type h 


Irs_cap(1,t) max num UAYs at LRS 1 in period t 
length(m) maximum periods of coverage of mission m 
val(m) value of mission m 
r(m) required periods of coverage for mission m 
k(s,Lt) indicator if Irs | used by schedule s in period t 
$OFFTEXT 
$OFFLISTING 
$INCLUDE uasData.gms 
PARAMETER 
k(s,1,0) 
$INCLUDE k.gms 
S$ONTEXT 
{DYNAMIC} SETS 
A(s,m,t) schedule s covers mission m in period t 
B(s,v,t) schedule s uses tier v UAV in period t 
Q(s,c,t) schedule s requires payload class c in period t 
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P(m,mp) mission prerequisite fixed mission m requires at least one 
mission mp in same period 


E(m,mp) exclusive missions m and mp cannot be accomplished in same 
period 


‘ 


{Notional data: will be replaced with INCLUDE statements and accompanying 
data files } 


S$OFFTEXT 
alias(m,mp); 
alias(v,vp); 
alias(t,tp); 

$INCLUDE a.gms 

$INCLUDE b.gms 
$INCLUDE q.gms 
$INCLUDE hp.gms 
$INCLUDE p.gms 
$INCLUDE mg.gms 
$INCLUDE mh.gms 
$INCLUDE gv.gms 
$ONLISTING 





SET sl(s,]); 

LOOP((s,)$(SUM(t,k(s,Lt))>eps), 
sl(s,l) = yes; 

); 


SCALARS errors,warnings ; 


errors=0 ; 
warnings=0; 
{Loops for data checking go here} 
IF(errors>0, 
PUT / ‘inconsistencies found in active and value data: ',errors:10:0/; 
68 


period 


PUT ‘(should I have shut you down here?' / ; 
2 


VARIABLE 
OBJECTIVE 
BINARY VARIABLES 
X(s) Schedule s flown by some UAV 
Y(m,v,t) Mission m covered in period t 
W(im,v,g,t) Mission m handled by GCS g in period t 
INTEGER VARIABLES 
D(m) Total number of dwell periods on mission m 
MODS (c,v,l) 
EQUATIONS 


eqnObj {RO} Objective function measures total value of mission 
accomplishment 


UAVLimit(l,v,t) {R1} limit on UAVs flying by tier and period 
ModuleLimit(c,v,1,t) {R2} limit on payload modules flying by payload and 


TotalMods(c,v) {R2.5} 
GCSLimit(g,t,h) {R3} limit on GCS use by gcs period and type 
LRSLimit(1,t) {R4} limit on LRS use by LRS and period 
CoverageLimit(m,v,t) {R5} mission and period coverage control 
DwellLimit(m) {R6} limit on total dwell time by mission 
Prerequisite(m,v,t) | {R7} prerequisite controls by period 
Mission_Control(m,v,t) 

* Exclusive(m,mp,t) | {R8} mutually exclusive missions by period 


RequiredDwell(m) {R9} required missions must be completely covered 
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eqnObj.. {RO} 

OBJECTIVE =e= sum(m, val(m)*D(m)) 
UAVLimit(L,v,t).. {Rl} 

SUM(s$(B(s,v,t) and sl(s,l)), X(s)) =l= num_uavs(l,v,t) 
ModuleLimit(c,v,],t).. {R2} 

SUM(s§(sl(s,l) and B(s,v,t) and Q(s,c,t)), X(s)) =l= MODS(c,v,1) 
TotalMods(c,v).. 

SUM(, MODS (c,v,1)) =l= num_mods(c,v) 
GCSLimit(g,t,h).. {R3} 


SUM((m,v)$(MG(m,g) and MH(m,h) and GV(g,v)), W(m,v,g,t)) == 
gcs_cap(g,t,h) 


LRSLimit(1,t).. {R4} 

SUM(s$(k(s,1,0>eps), X(s)) =l= Irs_cap(,t) 
CoverageLimit(m,v,t).. {R5} 

Y(m,v,t) =l= SUM(s$(A(s,m,t) and B(s,v,t)), X(s)) 
DwellLimit(m).. {R6} 

D@m) =l= SUM((v,t), Y(m,v,t)) 
Prerequisite(m,v,t)$HP(m).. {R7} 

Y(m,v,t) == SUM((mp,vp)$P(m,mp), Y(mp,vp,t)) 
Mission_Control(m,v,t)$(not HP(m)).. 

Y(m,v,t) == SUM(g$(MG(m,g) and GV(g,v)), W(m,v,g,t)) 


* Exclusive(m,mp,t)$E(m,mp).. {R8} 
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* “Y(m,t) + Y(mp,t) =Il= 1 


ke 


RequiredDwell(m)$(r(m)>eps).. {R9} 
D(m) =g= r(m) 

LOOP(m, 
D.up(m) = length(m); {R10} 

iy 


MODEL UAS_PLANNER / 
eqnObj 
UAVLimit 
ModuleLimit 
TotalMods 
GCSLimit 
LRSLimit 
CoverageLimit 
DwellLimit 
Prerequisite 

* Exclusive 
Mission_Control 
RequiredDwell 

/ 


’ 


IF(errors=0, 

PUT 'solve UAS_-PLANNER'/ ; 
* UAS_ PLANNER.reslim=60. ; 
* UAS_ PLANNER. optcr=0.05 ; 


PUT ' reslim=',UAS_PLANNER.reslim:5:0,' seconds, 
optcr='",UAS_PLANNER.optcr:5:2 ; 


PUT ' optfile='",;UAS_PLANNER.optfile:2:0 / ; 


a 


SOLVE UAS_PLANNER USING MIP MAXIMIZING OBJECTIVE; 


PUT "PLANNER solver and model status 
"UAS_PLANNER.solvestat:4:0,UAS_PLANNER.modelstat:4:0 / ; 


IF( UAS_PLANNER.modelstat<>1 and UAS_PLANNER.modelstat<>8, 
errors=1 ; 
ELSE 
PUT / ‘Objective: ',;OBJECTIVE.1:10:4 /; 
PUT / ‘Module Deployment / ; 
LOOP((c,v,1)$(MODS.L(c,v,1l)>eps), 
PUT ' ‘/c.tl:10,' ',v.tl:8,' ',1.t1:10,': ,MODS.L(c,v,1):4:0 / ; 
)3 
PUT /'UAV flight schedules..." // ; 
LOOP(s§$(X.1(s)>eps), 
PUT' ‘\s.tl:10,' '; 
LOOP(v$SUM(t$B(s,v,t), 1), 
PUT' ‘v.th:8 ; 
); 
PUT' '"/.X.1(s):4:1/; 
LOOP(m$(SUM(t$A(s,m,t), 1)>0), 
PUT' '‘,mtl:10; 
LOOP(t$A(s,m,t), 
PUT ''/t.tl:5 ; 
)3 
PUT/; 
) 
); 
PUT /'GCS Usage'/ ; 
LOOP(g$(SUM((m,v,t)$(MG(m,g) and GV(g,v)),W.1(m,v,g,t))>eps), 
PUT' '7,g.tl:8/; 
LOOP((m,v)$(MG(m,g) and GV(g,v) and (SUM(t,W.1(m,v,g,t))>eps)), 
PUT ‘"»,;m.th:8,' '\(SUM(t,W.1(m,v,g,t))) ; 
LOOP(t$(W.1(m,v,g,t)>eps), 


AZ 


PUT ''/t.tl:8 ; 
)3 
PUT/; 
); 
); 
PUT /'MISSION COVERAGE '/ ; 
PUT' Mission Dwell Time'/; 
LOOP(m$(D.1(m)>0), 
PUT! ';m.tl:8,' ',D.1Gm):5:1/; 
LOOP (t$(SUM(v, Y.1(m,v,t))>eps), 
PUT'  '/t.th:4,' "((SUM(v, Y.1(m,v,t))):4:1 / ; 
); 
)3 
PUTCLOSE UASP_LOG; 


PUT MISSIONS_CSV ; 
PUT ',,, Total, Total, Value,Total' /; 
PUT ',Start,End,Time Periods,Time Periods,Per Time, Value' /; 


PUT 
"Mission, Time, Time,Requested, Assigned,Period,Achieved,',OBJECTIVE.L /; 


LOOP(m, 
PUT m.1l,’,,,,',D.1(m),’,',val(m),',',(D.1@m)*val(m)) /; 
); 
PUTCLOSE MISSIONS_CSV; 
{Rest of report writer goes here } 
); 
); 


IF(errors>0, 


PUT ' <*** errors aborting this run:’,errors:10:0 / ; 
Vs 


IF(warnings>0, 
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PUT /' <+++ warnings generated by this run: ',;warnings:10:0 / ; 
PUT ' please find out what these indicate' / ; 

); 

PUTCLOSE UASP_log ; 


IF(errors=0, 
{Build CSV Output Files} 
PUT UASP_STA; 
PUT ‘Optimal solution found. Successful completion. /; 


PUTCLOSE UASP_STA; 
) 
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